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Abstract 


A  theory  is  presented  which  can  be  used  to  investigate  the 
dynamic  response  of  arches  subjected  to  cyclic  support  movement.  The 
theory  assumes  that  the  arch  is  circular  having  a  constant  mass  per 
unit  arch  length,  the  cross-section  is  constant,  damping  is  small  and 
the  center  line  of  the  arch  is  inextensible.  The  effects  of  rotatory 
inertia  and  shear  deflections  are  assumed  small,  compared  with  the 
effects  of  flexure,  and  are  neglected.  Only  deflections  in  the  plane 
of  the  initial  curvature  of  the  arch  are  considered.  A  closed  form 
solution  is  obtained  for  free  vibrations  and  an  eight- term  series  in 
terms  of  the  free  modes  of  vibration  is  used  to  evaluate  the  steady 
state  solution  for  forcing  frequencies  up  to  the  fourth  free  flexural 
mode.  Experimental  evaluation  of  the  theory  is  presented  on  the 
basis  of  radial  displacements  for  arches  having  pinned-end  boundary 
condi  tions . 

Results  are  presented  for  arches  having  radius  to  thickness 
ratio  from  121.2  to  179.4  and  half-opening  angles  of  97.2  to  119.2  degrees 
for  forcing  frequencies  up  to  the  fourth  flexural  mode  of  free  vibration. 
For  some  arches  tested,  out-of-plane  displacements  were  noted  to  occur 
at  frequencies  having  values  approximately  sixteen  times  the  value  of 
the  frequency  for  the  first  free  flexural  mode.  A  discussion  of  the 
parameters  that  influence  the  agreement  between  the  theoretical  and  the 
experimental  results  is  included  in  the  presentation  of  experimental  data. 
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Notation 


All  symbols  are  defined  when  used,  however,  a  summary  of  the 
more  commonly  used  symbols  are  listed  below, 
constant 

time  dependent  function  for  steady  state  solution 
constant  in  the  steady  state  solution 
E  Young's  Modulus 

G  amplitude  of  support  movement 

g(t)  function  describing  support  movement 

H  arch  thickness 


cij 

Dj(t) 

d: 

J 


I 

k 

m 

PR(4>,t) 

PTU,t) 

R 

R 

rn 

t 

u(cf>,t) 

up(<M) 


second  moment  of  area  of  the  cross-section  about  the  centroidal 
axi  s 

radius  of  gyration 

mass  per  unit  length  of  arch 

external  radial  load  per  unit  arch  length 

external  tangential  load  per  unit  arch  length 

radius  of  curvature 

Rayleigh  quotient 

roots  to  an  auxiliary  equation 

time  coordinate 

radial  displacement 

radial  steady  state  displacement 


x 


U  (<f>)  n-th  radial  mode  shape 

U^(<|>)  absolute  steady  state  radial  displacement 

U^/G  radial  magnification  ratio 

v(4>, t)  tangential  displacement 

v  (4>,t)  tangential  steady  state  displacement 

V  (d>)  n-th  tangential  mode  shape 

V^(<J>)  absolute  steady  state  tangential  displacement 

Va/G  tangential  magnification  ratio 

w(<j>,t)  represents  either  v(<j),t)  or  u(4>,t) 

W  (<f>)  represents  either  V  (<j>)  or  U  (<j>) 

(  )*,(  )**,  etc.  means  »  etc. 

(  ) ,  (  )  means  ^ 


Greek  Symbols 

a  one-half  the  arch  opening  angle 

angular  coordinate 


A 


n 


00 

n 


n 


K 

n(  t) 


0)  /  p  T  ** 

n-th  dimensionless  frequency  (A  =  -5-/  —  ) 

n  ^  m 

n-th  natural  frequency 
frequency  of  support  movement 
curvature  (1/R) 
time  dependent  function 
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Object  and  Scope 


The  purpose  of  this  study  was  to: 

(a)  develop  and  present  a  theoretical  method  for  evaluating  the 
dynamic  response  of  circular  arches  subjected  to  cyclic  support  move¬ 
ment,  and 

(b)  experimentally  evaluate  the  theory. 

To  simplify  the  analysis,  a  circular  arch  with  a  constant 
cross-section  and  constant  mass  per  unit  length  was  chosen.  Only 
flexural  deflections  in  the  plane  of  the  initial  curvature  of  the  arch 
were  studied.  In  the  theoretical  analysis,  a  closed  form  solution  for 
the  governing  differential  equations  was  obtained  for  arches  having 
various  types  of  boundary  conditions.  The  experimental  investigation 
was  restricted  to  arches  having  pinned  ends. 
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Chapter  I 

Introduction  and  Literature  Review 


1 .1  Introduction 

This  thesis  consists  of  a  review,  in  Chapter  I,  of  pertinent 
previous  investigations  dealing  with  deep  arches  i.e.  opening  angle 
greater  than  80  degrees,  that  have  been  reported  in  the  literature.  On 
the  basis  of  this  review,  the  theoretical  analysis  is  summarized  in 
Chapter  II  and,  where  necessary,  expanded  to  meet  the  objectives  of 
this  investigation.  Chapter  III  contains  a  discussion  of  the  experi¬ 
mental  technique  and  a  description,  of  the  experimental  apparatus  used. 

A  comparison  between  theoretical  and  experimental  results  for  four 
arches,  having  radius  to  thickness  ratios  in  the  range  121.2  £  q- £  179.4 
and  one-half  opening  angles  in  range  97.2  £a  £  119.2  degrees,  for  pinned- 
end  boundary  conditions  is  also  presented  in  Chapter  III.  Chapter  IV 
contains  a  summary  of  the  more  important  observations  arising  from  the 
investigations  and  a  list  of  some  associated  problems  suggested  for 
further  study. 

The  coordinate  system  and  components  of  displacement  of  a 
point  on  the  middle  surface  of  the  arch  are  shown  in  Fig.  1.1-1. 

The  notation  used  is: 
b  arch  width 

E  Young 's  Modul  us 

H  arch  thickness 


2 


<J>=0 


(})=a 


FIG.  1.1-1  COORDINATE  SYSTEM  AND  COMPENENTS  OF  DISPLACEMENT 
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m 

R 

t 

u(<f>,t) 

v(  4> ,  t) 


a 

4> 


second  moment  of  area  of  the  cross-section  about  the 

centroidal  axis 

mass  per  unit  length  of  arch 

radius  of  curvature 

time  coordinate 

radial  displacement  measured  from  the  undeformed  shape  of 
the  arch 

tangential  displacement  measured  from  the  undeformed  shape 
of  the  arch 

one-half  the  arch  opening  angle 
angular  coordinate 
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1 .2  Literature  Review 
Lamb  ( 1 ,  1888) 

Considering  an  element  of  an  arch  and  using  a  combination  of 
potential  energy  techniques  and  calculus  of  variations,  Lamb  obtained 
the  following  equation: 


2 

mR  3<J> 


4  2  4 

3  v  ,  3  v,  ,  ,  3  v 

4  +  2>  +  (  2 

3<p  3<J>  34>  3t 


1 

m 


3P 


R 


3<j> 


1.2-1 


where  P^.(cj),t)  and  P^(4>,t)  are  external  loads  acting  in  the  tangential 
and  radial  directions  respectively.  In  deriving  eq.  1.2-1  Lamb  assumed 
the  mid-surface  of  the  arch  to  be  inextensible  which  relates  u  and  v 
through  the  equation 


u 


1  .2-2 


For  the  case  of  free  vibration  (PR  s  P^  e  0),  Lamb  assumed 
the  tangential  displacement  as 

v  =  l  VU)  e1  “nt  .  1.2-3 

n=l 

Substituting  eq .  1.2-3  into  1.2-1,  yields 


v!1  +  2vJ,V  +  v“(l-X$)  +  xl  Vn  =  0  1.2-4 

n  n  n  n  n  n 

,  2  _  2  mR4 

An  °°n  El 


where 


1  .2-5 
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and  the  notation  (  )*,  (  )**,  etc.  denotes  differentiation  with  respect 
to  the  angular  coordinate  <j>.  When  the  assumed  solution 
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l 

j=l 


was  substituted  into  eq.  1.2-4,  Lamb  obtained 


rn  +  2rn  +  rn  (1_An)  +  =  0  • 

n  n  n  n  n 


1.2-6 


Lamb  did  not  obtain  a  solution  to  eq.  1.2-6  for  the  general  case,  but 

he  did  obtain  some  results  for  a  shallow  arch  by  noting  that  one  root 

2 

of  the  cubic  in  rn  is  nearly  equal  to  unity  and  approximating  the  re- 

2 

maining  roots  by  rn  +  1  =  +  A.  With  this  approximation,  the  mode  shapes 
and  the  natural  frequencies  for  the  lower  modes  were  obtained. 

By  neglecting  the  inertia  terms  and  external  loads  in  eq. 

1.2-1  Lamb  also  solved  for  the  deflection  of  arches  under  various 
static  loads  acting  at  the  boundaries,  <|)  =  +  a,  of  the  arch. 

Den  Hartog  (2,  1928) 

By  use  of  the  Rayleigh  quotient.  Den  Hartog  calculated  the 
lowest  natural  frequency  of  circular  arches  subtending  any  angle  up  to 
360  degrees,  for  pinned  and  clamped  ends.  For  the  tangential  displace¬ 
ment  he  assumed  a  function 
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v  =  —  (cos  ^  +  1 )  sin  cot 

it  a 

and  by  satisfying  inextensibility,  the  radial  displacement  was 

u  =  A  sin  ^  sin  oot 

a 

for  pinned-end  boundary  conditions.  Using  the  above  displacements, 
the  maximum  strain  and  kinetic  energies  were  equated  and  the  natural 
frequency  oo  was  obtained.  Den  Hartog  stated  that  for  30°  <  a  <  135°, 
displacements  of  the  form 

and  v  *  Aa  -  (^)2  +  —  (^)4  -  j-  (^)6]  sin  cot 

yield  smaller  and  therefore  better  results  for  go.  Upon  checking  the 
frequency  equation  Den  Hartog  obtained  by  using  the  polynomial  displace¬ 
ments,  an  error  was  noted,  and  in  fact  the  polynomial  solution  yields 
higher  and  therefore  poorer  values  for  the  natural  frequency  than  does 
the  assumed  sinusoidal  displacement  for  all  values  of  a.  The  equation 
Den  Hartog  obtained  was 


A2  = 


1 

T 

a 


a4 -78. 5a2  +1585 

1+0.094  a2 
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whereas  the  equation  should  have  read 

,2  _  1  a4-79.20a2+1584 
A  4  T~  • 

a  1+0.07769  a 

Using  an  analysis  similar  to  that  above,  Den  Hartog  also 
evaluated  natural  frequencies  for  fixed-end  boundary  conditions.  Den 
Hartog  investigated  extensional  arches,  noting  that  for  certain  R/k 
(k  =  radius  of  gyration  of  the  cross-section)  and  a  values,  the  ex¬ 
tensional  mode  had  a  lower  natural  frequency  than  the  inextensional 
mode . 

Love  (3,  1944) 

Love  used  the  classical  approach  to  obtain  the  governing 
differential  equation  for  a  thin  arch  by  applying  conditions  of  dynamic 
equilibrium  to  an  el  enent  cut  from  the  arch  as  shown  in  Fig.  1.2-1. 


Q 


FIG.  1.2-1  ELEMENT  CUT  FROM  ARCH 
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Since  Love  investigated  free  vibration,  the  external  loads  PD  and  PT 

K 

are  zero  for  this  case.  The  three  equations  obtained  were: 


|3-  +  N  -  mR  Ajj.  *  0 

34  3? 


1  .2-7 


3N  n  mD  3  v  _  n 
-  Q  -  mR  — n-  -  0 

34  3? 


1  .2-8 


and 


1  AM  =  Q 

R  3<j)  ^ 


1,2-9 


where  rotatory  inertia  and  shear  deflection  effects  were  assumed  to 
be  negligible  compared  with  the  effects  of  flexure.  The  bending  moment 
M  and  normal  force  N,  expressed  in  terms  of  the  displacements,  are 


and 


m  -  -  4-  (^4+  u) 


7 


w 


1  .2-10 


,  3v  v 
34>  "  U 


1 .2-11 


For  inextensibility  of  the  center  line,  the  tangential  and  radial  dis¬ 
placements  are  related  by 


3v_  _ 

3<J>  ' 


u 


1  .2-12 


Combining  eqs .  1.2-7  through  to  1.2-12,  Love  obtained  the  equation 


8 


El  ,36w 

~T  6 

mR  3(f) 


34w 

2  2 

3(f)  3t 


32w  _ 


=  0 


1  .2-13 


where  w  represents  either  the  tangential  displacement  v  or  the  radial 
displacement  u.  The  solution  to  eq.  1.2-13  was  of  the  form 


w((j),t)  =  l  W  ( 4>)  cos  (wnt  +  0)  1.2-14 

n=l 


which  upon  substituting  into  eq.  1.2-13,  resulted  in  an  equation,  with 
one  independent  parameter  <J>,  of  the  form 

W ll  +  2W*V  +  wJ/O-A2)  +  A2  =  0  1.2-15 

n  n  n  n'  n  n 

2 

where  An  is  given  by  eq.  1.2-5.  A  comparison  of  eq.  1.2-15,  as  obtained 
by  Love,  with  eq.  1.2-4,  as  obtained  by  Lamb,  reveals  that  they  are 
identical.  A  solution  for  W  (<f>)  was 

3 

Wn<f>  *  ^  (Ank  cos  rnk^  +  Bnk  sin  rnk*> 


where  r  is  obtained  from 
n 

r2  (r2-!)2  -  (r2+l)  A2  =  0  1  .2-16 

n  n  n  n 

which  is  identical  to  eq .  1.2-6  as  obtained  by  Lamb.  For  a  complete 
ring  rR  must  be  an  integer,  and  there  are  vibrations  with  r^  wavelengths 


9 


to  the  circumference,  rn  being  any  integer  greater  than  unity.  The 
natural  frequencies  of  vibration  for  a  complete  ring  are  given  by 


co_  = 


rV-l)2 

El  nirn  ' 

— W  — - 

mR  r^  +  1 


1  .2-17 


Hoppe  (4)  is  credited  for  obtaining  eq.  1.2-17  in  1871. 

Love  suggests  a  method  for  obtaining  values  of  r^  which  deter¬ 
mine  the  natural  frequencies  of  vibration  for  complete  rings,  however, 
no  values  of  rn  for  arches  with  specific  opening  angles  and  boundary 
conditions  were  recorded. 


Bolotin  (5,  1952;  6,  1953;  25,  1956) 

Using  the  basic  equations  developed  by  Lamb,  Bolotin  solved 
the  stability  problem  of  arches  under  parametric  excitation  at  the  crown 
(4)  =  0)  of  the  arch.  Lamb's  equations  were  modified  to  account  for  the 
dynamic  load  acting  at  the  crown.  Bolotin  investigated  this  problem 
experimentally  obtaining  good  agreement  between  theory  and  experiment, 
however,  the  experimental  details  were  not  given.  Another  aspect  of 
Bolotin's  work  (25)  involved  investigating  the  dynamic  response  of 
arches  and  rings  subjected  to  pulsating  pressure  loads. 

Morley  (7,  1957) 


Morley  calculated  the  natural  frequencies  and  normalized  mode 
shapes  for  the  first  20  free  modes  of  vibration  for  a  thin  elastic  arch 
subtending  an  angle  of  360  degrees  having  free-end  boundary  conditions. 
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Morley's  approach  to  the  cut-ring  problem  was  identical  to  the  method 
used  by  Love  to  solve  for  the  modal  characteristics  of  a  complete  ring. 
However,  since  the  ring  was  cut,  the  roots  of  eq.  1.2-16  are  no  longer 
integers  but  are  determined  by  satisfying  the  free-end  boundary  condi¬ 
tions.  In  satisfying  the  boundary  conditions,  Morley  solved  for  the 
eigenvalues  of  a  matrix  from  which  it  was  possible  to  calculate  the 
natural  frequencies.  The  eigenvectors  associated  with  these  eigen¬ 
values  determined  the  mode  shapes.  Results  were  not  presented  for  any 
other  boundary  conditions  or  arch  opening  angles. 

Archer  (8,  1959) 

In-plane,  inextensional  vibrations  of  an  arch  were  studied  by 
Archer,  using  the  basic  equations  as  derived  by  Love,  with  the  addition 
of  terms  to  represent  damping  effects.  This  equation  can  be  represented 
as 


El  ,  36v  ?  34v 

— t  \ — g-  +  t  — t-  + 

mR  3<j>  34> 


32v.x 

3<T 


34V 

3(j)23t2 


32v 


3t 


+  £( 

m  v 


33v 


3t3(f> 


I 


3Vx  _ 

"  3t' 


0  1.2-18 


where  c  represents  the  viscous  damping  coefficient.  Archer  solved  eq. 
1  .2-18  by  using  a  product  solution  of  the  form 


v(4>,t)  =  VU)  T(t) 


1  .2-19 


which  separated  eq.  1.2-18  into 


. 
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mK 


1  .2-20 


and 


1 .2-21 


Love's  approach  was  used  to  reduce  eq.  1.2-20  to  eq.  1.2-17.  In  a 
similar  approach  to  that  of  Morley,  Archer  solved  for  the  roots  of  eq. 

1 .2-17. 

Archer  extended  the  solution  of  eq.  1.2-20  to  find  the  transient 
response  of  a  semi-circular  arch  under  a  small  step  motion  of  one  support 
relative  to  the  other.  The  particular  problem  was  defined  by  assuming 
one  support,  say  =  a,  was  initially  displaced  tangentially  by  a  small 
amount  R  ,  where  Aip  defines  the  change  in  slope  between  the  deformed 
and  undeformed  state  at  the  support.  The  end  was  now  allowed  to  return 
to  its  natural  position  according  to  the  motion 


v(a,t)  =  R  AiJj  e  ,  t  >  0  . 


1 .2-22 


Equation  1.2-22,  represents  a  time-dependent  boundary  condition.  In 
obtaining  a  solution,  Archer  employed  a  method  developed  by  Mindlin  and 
Goodman  (9)  for  solving  problems  with  time  dependent  boundary  conditions. 
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Eppink  and  Veletsos  (1  0,  1  959;  11  ,  1960;  12,  1960) 

Eppink  and  Veletsos  described  a  numerical  procedure  which  can 


be  used  for  computing  the  dynamic  response  of  circular  arches  subjected 
to  transient  forces  for  both  the  elastic  range  of  behavior  and  for  the 
range  approaching  failure.  The  problem  was  analyzed  approximately  by 
replacing  the  continuous  arch  by  a  framework  consisting  of  a  series  of 
rigid  bars  and  flexible  joints.  Any  external  loads  were  assumed  to  be 
concentrated  at  the  joints.  The  authors  evaluated  the  natural  frequencies 
and  corresponding  modes  for  a  two-hinged  and  for  a  fixed  arch  of  constant 
cross-section.  In  particular,  the  authors  considered  time-dependent 
pressure  as  the  external  load  acting  on  the  arch.  This  method  of  analysis 
permitted  a  calculation  of  the  critical  pressure  load  necessary  to 
buckl e  the  arch . 

This  work  was  extended  (11)  to  Include  the  behavior  of  arches 
in  the  inelastic  range  by  approximating  the  cross-sectional  area  of  the 
arch  by  two  flanges  connected  by  a  thin  rigid  web.  The  resistance  of 
each  flange  was  represented  by  a  bilinear  stress-strain  diagram.  The 
authors  obtained  results  for  a  uniform  all-around  pressure  with  a  time- 
dependent  variation  represented  by  a  triangle  with  an  initial  peak  and 
for  a  triangular  moving  pressure  pulse.  The  authors  also  investigated 
the  effects  of  load  parameters  such  as  duration  of  pressure  pulse, 
magnitude  of  peak  pressure  compared  to  the  critical  buckling  pressure 
and  the  velocity  of  propagation  of  the  pressure  front.  The  effects  of 
arch  parameters,  i.e.  geometric  and  physical  properties,  were  also  in¬ 
vestigated  . 
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The  previous  method  of  analysis  was  applied  to  a  two-hinged 
arch,  of  uniform  cross-section  and  mass  per  unit  length,  which  was 
subjected  to  a  pressure  uniformly  distributed  around  the  arch  (12). 

The  resulting  governing  equations  were  evaluated  on  a  computer  using 
a  step-by-step  method  of  integration.  These  results  and  the  results 
from  the  modal  method  were  compared  showing  good  agreement  between  the 
two  methods . 

Vo  1  terra  and  Morell  (13,  1960;  14,  1961;  15,  1961;  26,  1967) 

Vol terra  and  Morell  extended  the  Rayl eigh-Ritz  method,  used 
by  Den  Hartog,  to  find  the  lowest  natural  frequency  of  elastic  arches 
having  the  center-line  in  a  form  of  a  cyloid,  parabola,  catenary,  as 
well  as  a  circle  subjected  to  clamped-end  boundary  conditions  (13). 

This  investigation  was  extended  (14)  to  include  vibrations  outside  the 
plane  of  initial  curvature  of  the  arch.  A  correction  to  the  results  pre¬ 
sented  in  this  paper  (14)  was  given  in  (26).  For  hinged  boundary  con¬ 
ditions  (15)  the  lowest  natural  frequency  was  determined  for  arches 
vibrating  either  in  the  plane  or  out-of-the  plane  of  the  initial  curva¬ 
ture.  For  in-plane  vibrations,  they  obtained  the  natural  frequency  to 
be 
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1 .2-23 


where 
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F(a)  = 


“(RLsin2  (£t)2  (1  +  cos  ^)2])d<)> 
v  a  it  a  ' 


o 


For  a  particular  case  of  a  circle,  R  =  RQ,  eq.  1.2-23  reduces  to 


which  is  identical  to  the  result  obtained  by  Den  Hartog .  The  function 
F(a)  can  be  evaluated  for  a  cycloid,  catenary  and  parabola  by  sub- 


2  3 

stituting  for  R  the  values  R  =  RQ  cos  R  =  RQ/cos  <j>  and  R  =  Rq/cos  $ 


respectively  and  performing  the  necessary  integration.  These  authors 
presented  results  only  for  arches  with  maximum  opening  angles  of  up  to 
80  degrees . 

Lang,  T.E.  and  Reed,  R.E.  (16,  1  962) 


Lang  and  Reed  also  presented  an  approximate  method  for  deter¬ 


mining  in-plane  and  out-of -plane  natural  frequencies  and  mode  shapes 
for  arches  and  complete  circular  rings  with  non-uniform  stiffness  and 
mass  distribution.  The  ring  was  idealized  by  a  series  of  point  masses 
connected  by  massless  elastic  arch  segments.  The  transfer  matrix 
technique  was  used  to  obtain  a  solution.  These  results  were  compared 
to  the  closed  form  solution  obtained  by  Love  (3)  (Hoppe  (4),  1871)  for 
the  natural  frequencies  for  a  thin,  constant  thickness  complete  ring 
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with  resulting  differences  of  0,  1,  and  4.1  percent  for  the  second, 
fourth  and  sixth  mode  respectively,  for  in-plane  flexural  vibrations. 

For  out-of -plane  vibrations  the  differences  in  the  natural  frequencies 
were  0,  0,  1.5,  and  8.9  percent  for  the  second,  fourth,  sixth  and  eight 
mode  respectively.  In  all  cases  where  differences  occurred,  the  approxi¬ 
mate  technique  gave  lower  values  for  the  natural  frequencies  than  those 
obtained  from  the  closed  form  solution.  The  authors  attributed  the 
lower  values  to  the  fact  that  Hoppe  did  not  take  into  account  the  effects 
of  rotatory  inertia. 

Nelson  (17  ,  1962) 

A  modification  of  the  Rayleigh-Ritz  technique,  as  used  by  Den 
Hartog  (2),  was  made  by  Nelson  by  introducing  Lagrangian  multipliers  to 
take  into  account  the  boundary  conditions.  For  the  unsymmetrical  modes 
with  pinned-end  boundary  conditions,  the  radial  displacement  was  assumed 
to  be 


oo 


sin 


nTTd) 

2a 


which  was  used  together  with  the  Inextensibility  condition  given  by 
eq.  1.2-2  to  evaluate  the  energy  function  defined  as 


W  -  R  T 


where  W  is  the  strain  energy,  T  is  the "reduced  kinetic  energy'and  R  is 
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the  Rayleigh  quotient.  An  extremum  of  the  energy  function  was  then 
found  with  respect  to  the  A^'s  and  a  constant  B  which  was  obtained  from 
evaluating  the  tangential  displacement  V  by  using  eq.  1.2-2,  i.e. 

a 

V  =  |  U  d<()  +  B  . 

-a 

Introducing  the  Lagrangian  multiplier  y,  the  energy  function  was  written 
as 


T  I  [1 
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which  was  minimized.  For  a  one-term  approximation  to  the  resulting 
series  frequency  equation.  Nelson  expressed  the  nondimensional  frequency 
parameter  6 ,  as 


where  L  =  2Ra,  the  developed  length  of  the  arch.  The  roots  R-j ,  R^,  ... 

2 

obtained  from  the  minimizing  process  have  the  property  that  R-j  >_  oo-|  , 

2 

R^  ^  (^2  ,  ... 

Nelson  referred  to  the  phenomena  of  "frequency  cross-over" 
by  showing  that  the  arch  will  shift  from  a  one-node  to  a  three-node  mode 
of  vibration  by  increasing  the  arch  opening  angle.  Nelson  observed  that 
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this  occurred  at  an  opening  angle  of  about  570  degrees. 

Results  were  presented  for  arch  opening  angles  beyond  360 
degrees.  However,  in  the  theory,  it  was  assumed  that  the  arch  was  vi¬ 
brating  in  its  plane  of  initial  curvature.  To  satisfy  this  assumption 
it  is  physically  impossible  to  have  an  actual  arch  greater  than  360 
degrees.  Hence,  any  results  presented  by  Nelson  for  arch  opening  angles 
greater  than  360  degrees  have  no  physical  significance. 

Nelson  also  investigated  symmetrical  modes  and  extensional 
modes  by  using  the  above  method  of  analysis.  Graphical  results  were 
presented  indicating  for  what  arch  parameters  the  inextensional  theory 
yields  accurate  results  and  for  what  parameters  the  extensional  theory 
must  be  used. 


Lang  (18,  1962;  19,  1963) 


Except  for  Archer's  transient  solution  for  an  arch  when  one 
support  is  given  a  particular  displacement,  the  studies  conducted  thus 
far  consisted  of  finding  the  natural  frequencies  and  mode  shapes  for 
freely  vibrating  arches.  It  appears  that  Lang  was  the  first  to  propose 
a  closed  form  solution  for  arches  or  rings  which  are  excited  by  time 
dependent  external  forces.  A  theoretical  analysis  was  applied  to  both 
the  extensional  and  inextensional  case.  The  classical  approach  was 
used  to  derive  the  governing  differential  equations.  For  the  extensional 
case,  the  homogeneous  equation  obtained  was 


VVnI  +  (2+4?>VnV+n-^--4) 
n  (|)2  n  n  (£)  n 


+  *n  f1  -  TO)  Vn  ■  0 

{F 


1.2-24 
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where  k  is  the  radius  of  gyration  of  the  cross-section  and  the  terms 
2  2 

X  / (R/k)  take  into  account  the  effects  of  stretching  of  the  center  line. 

For  the  inextensional  case,  the  complete  governing  differential 
equation  obtained  was 


UL  I A  + 

4  L  r 

mR  d<$> 


2  4 

3Vi  .  9  v 

2J  2 
9<T  9(j)  9t' 
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m  ^  9<}> 


-  P-rU.t)] 


1.2-25 


where  and  Pj  are  external  loads  per  unit  arch  length  in  the  radial 
and  tangential  direction  respectively.  Also  for  the  inextensional  case, 
the  radial  displacement  u  and  the  tangential  displacement  v  are  related 
by  the  equation 

lx  =  u  .  1.2-26 

dtp 

Since  eq.  1.2-25  is  a  linear  equation,  the  steady  state  solution  Lang 
proposed  was  a  series  in  terms  of  the  free  modes  of  vibration,  i.e. 

oo 

vn(<M)  =  l  V.  {< (,)  n,(t)  1.2-27 

p  i=l  1  1 

where  is  the  i-th  mode  of  free  vibration  in  the  tangential  direction 
and  q.j(t)  is  a  function  of  time,  as  yet  unknown.  Substituting  eq. 

1.2-27  into  eq.  1.2-25,  multiplying  by  V .  dcp ,  and  integrating  between 

J 
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-a  and  a,  the  resulting  equation  is 

-  PT]d<t, 

1 .2-28 

2  2  III 

where  (“)  means  d  (  )/dt  ,  and  (  )  ,  (  )  ,  ...  means  d(  )/d<J>, 

2  2 

d  (  )/d<J>  ,  etc.  The  orthogonality  condition  for  the  free  modes  of 

vibration  is 


l  [  (n-  +  w-n-) 
i=l  1  1  1 


rot 
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(v’v'  +  V.V.)d*]  = 
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(vlvl  +  V.V.) d*  =  0  for  i  f  j 


and  for  i  =  j  the  modes  were  normalized  according  to 


_1_ 

2a 
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<v?v 5  +  viVd*  ■  1 


1.2-29 


Because  of  the  orthogonality  condition,  the  i-th  and  the  j-th  modes  un¬ 
coupled  in  eq.  1.2-28  and  using  eq.  1.2-29,  the  resulting  equation  for 
n j ( t )  was  found  to  be 


n.  +  u.n.  ■  Dj(t) 


1.2-30 


a 


where 


V,  8P 


R 


j  m  (^r1  pT) 
D.(t)  =  -  - 1 - 


Lang  obtained  the  steady  state  solution  in  integral  form  but  did  not 


20 


evaluate  it  for  any  particular  type  of  external  loads. 

Experimental  evaluation  of  the  homogeneous  solution  to  eq. 
1.2-25  was  obtained  showing  good  agreement  between  theory  and  experiment 
for  the  natural  frequencies  and  corresponding  mode  shapes  for  the  first 
three  symmetric,  in-plane  modes  for  a  complete  ring  with  a  point  clamp¬ 
ing,  i.e.  the  ring  was  essentially  an  arch  with  fixed  ends  and  a  =  180 
degrees.  No  experimental  verification  was  presented  for  other  arches. 

In  a  later  report,  Lang  (19)  tabulated  the  theoretical  di¬ 
mensionless  frequency  parameters  as  given  by 

.2  2  mR4 

Xn  =  “n  ET 


and  the  corresponding  mode  shapes  for  the  first  five  in-plane  modes 
for  arches  with  opening  angles  of  180  degrees  for  various  types  of 
boundary,  conditions. 

Takahashi  (20,  1963) 

Lagrange's  equation  was  used  by  Takahashi  to  obtain  the 
governing  differential  equation  for  the  vibration  of  a  circular  arch 
in  its  initial  plane  of  curvature.  This  approach  yielded 
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1.2-31 


-a 


where  sinusoidal  time  functions  were  assumed  and  L  is  the  Lagrangian. 

2 

The  value  of  An  is  given  by  eq.  1.2-5.  The  method  of  calculus  of  vari¬ 
ation  was  applied  to  eq.  1.2-31  and  upon  rearranging  the  resulting 
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equations  the  following  equation  was  obtained 
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1.2-32 


This  equation  is  identical  to  the  equation  Lang  obtained  in  his  analysis 
(eq.  1.2-24).  If  arch  parameters  are  chosen  such  that 


i.e.  stretching  of  the  mid-surface  is  negligible,  eq.  1.2-32  reduces 

to  the  equation  Lamb  obtained  (eq.  1.2-4).  Takahashi's  solution  in- 

o 

dicated  that  the  slenderness  ratio,  ( R/k)  ,  influences  the  resonant 
frequencies  for  symmetric  modes,  however,  this  ratio  has  little  effect 
on  the  frequencies  for  unsymmetrical  modes.  Takahashi  also  observed 
the  "frequency  cross-over"  reported  by  Nelson  (17). 


Nelson  (21  ,  1964) 

Extending  his  earlier  work  (17),  Nelson  investigated  the 
frequency  cross-over  between  extensional  and  i nextensional  modes  for 
pinned-end  circular  arches.  Results  were  presented  showing  the  effects 
of  the  parameter,  g  =  (L/k-rr)  ,  on  the  mode  of  vibration,  where  L  is 
the  developed  length  of  the  arch.  These  results  indicated  that  for 
certain  arch  opening  angles  and  certain  values  of  g,  the  arch  vibrated 
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in  the  extensions!  mode  rather  than  in  the  inextensional  mode.  It  was 
shown  that  there  was  a  unique  value  of  the  frequency  where  this  cross¬ 
over  would  occur  which  Nelson  called  the  cross-over  point. 

Another  aspect  of  Nelson's  investigation  involved  finding 
over  what  range  of  a(a  =  one-half  the  arch  opening  angle)  and  for  what 
values  of  g  the  inextensional  theory  can  be  used  to  obtain  sufficiently 
accurate  results. 

Antman  and  Warner  (22,  1965) 

Starting  with  the  basic  equations  as  established  by  Love, 
Antman  and  Warner  investigated  the  stability  of  vibrating  arches.  The 
governing  equations  were  written  as 


1.2-33 


and 


-  k  N  -  Qq  +  m  ^  =  0 


1.2-34 


where  Mq,  Nq  and  Qq  are  the  externally  applied  moment,  tangential  load 
and  shear  load  respectively  per  unit  length  of  arch,  k  is  the  deformed 
curvature  and  ds  is  an  element  of  arch  length.  Antman  and  Warner  argued 
that  if  the  moment  M,  normal  force  N  and  curvature  k  were  expressed  in 
terms  of  the  displacements  u  and  v,  the  term  kN  in  eq.  1.2-34  gave  rise 
to  a  nonlinear  coupling  in  displacements.  However,  if  N  was  regarded 
as  a  known  function  of  s  and  t,  then  kN  contributed  a  linear  term  to 
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eq.  1.2-34.  Thus,  these  authors  obtained  an  independent  solution  for 
N  from  the  linear  equations  and  substituted  this  value  of  N  into  eq. 
1.2-34  which  was  then  examined  for  stability.  For  a  circular  arch  of 
radius  R,  the  resulting  equation  was 
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where  N  was  obtained  from 
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for  the  extensional  case  and  from 


a2N  8Nn  a2., 
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1.2-37 


for  the  inextensional  case.  Although  N  could  be  solved  from  eq.  1.2-36 
or  eq.  1.2-37,  the  authors  obtained  N  by  linearizing  the  governing 
differential  equations  by  assuming  a  constant  curvature  before  and  after 
deformation  and  then  solving  for  a  steady  state  value  of  N  using  Lang's 
(18)  approach.  Denoting  the  steady  state  value  of  N  by 


N  =  N*(<j))  cos  fit 


and  substituting  this  value  of  N  into  eq .  1.2-35  the  following  equation 
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was  obtained 


1.2-38 


where 


f(4>)  ■  +  QgU)  -  l 
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and 


Q0  =  Qg(4>)  cos  fit  ,  =  Mg(4>)  cos  fit  . 


The  problem  was  now  reduced  to  finding  under  what  conditions  the  solu¬ 
tion  u  to  eq.  1.2-38  would  remain  bounded  for  all  time,  t  >  0.  Since 
eq.  1.2-38  was  not  separable,  a  Kantorovich's  modification  of  Galerkin's 
method  to  reduce  eq.  1.2-38  to  a  system  of  coupled  Mathieu  equations 
was  employed. 


Antman  and  Warner  analyzed  the  stability  of  a  semi-circular 


arch  under  radial  forced  end  vibrations  as  shown  in  Fig.  1.2-2.  From 
this  investigation  it  was  concluded  that  "secondary  resonant  frequencies" 
were  given  by 
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where  j  is  any  integer,  i  may  equal  k  and  $k  is  defined  as 

Bk  =  ^7  (k+1)2  [(k+1>2  -  1]  • 
mR 


FIG.  1.2-2  ARCH  MODEL  USED  BY  ANTMAN  AND  WARNER 


Popescu  (27,  1970) 

To  investigate  the  dynamic  stability  of  arches,  Popescu  used 
an  energy  method  similar  to  that  used  by  Bolotin  to  obtain  the  govern¬ 
ing  equation.  The  arches  were  assumed  to  be  subjected  to  a  distributed 
load  q(t)  given  by  the  equation 

q ( t )  =  Qq  +  Q  cos  fit 
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where  Qq  is  the  static  radial  load,  Q  is  the  amplitude  of  the  periodically 
variable  load  and  ft  is  the  pulsating  frequency.  The  governing  equation 
was  reduced  to  a  Mathieu  equation  which  was  then  analyzed  for  stability. 

In  obtaining  the  expression  for  the  kinetic  energy,  Popescu  did  not 
account  for  the  tangential  displacement  of  the  arch.  This  resulted 
in  a  higher  bound  for  the  frequency  of  free  vibration  than  that  ob¬ 
tained  by  Den  Hartog,  Nelson,  Eppink  and  Veletsos,  and  Bolotin  who  con¬ 
sidered  the  tangential  displacement. 

1 .3  Closure  to  Chapter  I 

The  basic  equation  governing  free  vibrations  of  inextensible 
arches  had  been  obtained  by  Lamb  in  1888.  Since  that  time,  numerous 
approximate  techniques  were  employed  to  obtain  responses  of  deep  arches. 
For  example,  Den  Hartog  used  the  Rayleigh  quotient  to  obtain  the  charac¬ 
teristics  of  the  first  mode  for  both  the  extensional  and  inextensional 
case;  Bolotin  investigated  the  effects  of  parametric  excitation  and 
the  effects  of  pulsating  pressure  loads  on  arches  using  the  Galerkin 
method;  Eppink  and  Veletsos  obtained  the  dynamic  response  of  arches 
excited  by  dynamic  pressure  loads  by  replacing  the  continuous  arch  with 
a  framework  of  rigid  bars  and  flexible  joints;  Nelson  used  the  Rayleigh- 
Ritz  method  in  conjunction  with  Lagrangian  multipliers  to  calculate 
the  characteristics  of  the  free  modes  of  vibration;  Vol terra  and  Morel  1 
extended  the  Rayleigh-Ri tz  method  to  include  arches  with  center- lines 
other  than  circular;  and  Lang  and  Reed  used  the  transfer  matrix  technique 
to  solve  for  the  free  vibration  characteristics. 

Other  investigators  such  as  Archer,  Morley,  Love,  Hoppe  and 
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Takahashi  solved  the  governing  equation  for  arches  or  rings  for  free 
vibrations  in  closed  form.  Antman  and  Warner  considered  some  nonlinear 
terms  in  the  governing  equations  and  performed  a  stability  analysis  for 
arches.  It  appears  that  Lang  was  the  first  investigator  to  derive  the 
orthogonality  condition  for  the  free  modes  of  vibration  for  an  arch. 
Lang  proposed  a  steady  state  solution  as  a  series  in  terms  of  the  free 
modes  of  vibration.  This  steady  state  solution  was  stated  in  integral 
form  and  was  not  evaluated  for  any  type  of  loading. 

Very  little  experimental  data  is  available  to  evaluate  the 
above  approximate  and  closed  form  solutions.  Lang  conducted  experi¬ 
mental  tests  on  complete  rings  only,  obtaining  good  agreement  between 
theory  and  experiment  for  the  natural  frequencies  and  mode  shapes. 
However,  no  experimental  evaluation  was  provided  for  the  steady  state 
solution. 

On  the  basis  of  the  literature  review,  it  is  evident  that  an 
extension  to  the  presently  existing  theory  is  necessary  in  order  to 
provide  a  solution  for  the  dynamic  response  of  arches  subjected  to 
cyclic  support  movement.  The  theoretical  analysis  to  the  above  problem 
is  presented  in  Chapter  II. 
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Chapter  II 
Theory 

2.1  Derivation  of  the  Governing  Differential  Equations 

The  governing  differential  equation  used  in  this  analysis 
will  be  derived  by  considering  the  equilibrium  of  an  arch  element. 
Lang's  (18)  approach  for  obtaining  the  steady  state  solution  will  be 
modified  to  account  for  the  inertial  loads  caused  by  support  movement. 

A  steady  state  solution  for  symmetrical  and  unsymmetrical  support  move¬ 
ment  will  be  proposed.  The  coordinate  system  used  is  shown  in  Fig. 
1.1-1. 

The  assumptions  used  in  the  theory  are: 

1.  The  material  is  perfectly  elastic  and  obeys  Hooke's  Law. 

2.  The  displacements  u(<j>,t)  and  v(cj),t)  measure  deformation 
only,  i.e.  rigid  body  motion  is  not  accounted  for  by  u  and  v. 

3.  The  mass  per  unit  arch  length  m,  thickness  H,  and  width  b 
are  constant  along  the  arch  circumference. 

4.  Displacements  out-of-the  plane  of  initial  arch  curvature 
are  not  considered. 

5.  Although  the  existence  of  damping  is  acknowledged,  it  is 
assumed  that  its  effects  on  the  steady  state  solution  are  negl igible. 

6.  The  arch  vibrates  with  no  extension  of  the  center  line 

i.e.  (R/k)2  »  . 
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7.  The  undeformed  shape  of  the  arch  is  circular  with  radius  R. 

8.  Radial  displacements  u  are  small  such  that  R  +  u  »  R. 

9.  Rotatory  inertia  and  shear  deformation  effects  are  small 
compared  with  the  effects  of  flexure. 

10.  Second  order  terms  are  negligible  when  compared  to  first 
order  terms. 

By  applying  dynamic  equilibrium  to  an  element  of  an  arch 
shown  in  Fig.  1.2-1,  the  three  governing  equations  are: 


2 

!§■  +  N  -  mR  t  R  P  =  0  , 

at2  R 


2.1-1 


It  -  Q  -  mR  +  R  PT  =  0  , 

^  at2  T 


2.1-2 


and 


la.n 

R  3<j>  w 


2.1-3 


where  (4> , t )  and  P^(<j),t)  represent  the  external  load  per  unit  arch 

length  acting  in  the  tangential  and  radial  directions  respectively. 
Combining  eq.  2.1-1,  2.1-2  and  2.1-3  together  with  eq.  1.2-10  and 
1.2-11  yields 


El  / 85u  ,  o  93u  ,  9Ux  83u  32v  =  1 

mR4  3(j)5  34>3  ^  3(})3t2  7?  m  8<f> 


2.1-4 


The  condition  of  inextensibility  as  given  by  eq.  1.2-12  separates  eq. 
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2.1-4  into 


EI  ,36v  ,,  34v  32V'  34v 

4  ^  6  +  ^  4+  2'  +  22 

mR  3<J>°  9<p  9(f)  94>  at 


32V  .  1  ,3PR 

17  ’ m  34 


-  PT> 


2.1-5 


and 


6  4  2  4 

EI  , 3%  2  £u  +  £u w  d\ 

T  \  g  r  c  4t  2'  22 

mR  3(f)0  9cf)  9<JT  9<f>  9t 


32u  =  1  / 
m  ' 


2 

9  K 


R 


9t 


7 


3<p 


9R1 

9<J) 


-)  2.1-6 


It  will  be  shown  that  eq.  2.1-6  reduces  to  the  governing 
equation  of  a  straight  beam  in  the  limit  as  R  •*  «>.  Rewriting  eq.  2.1-6 
using  an  element  of  arch  length  ds  as  the  independent  variable  where 
Rdcj)  =  ds ,  yields 


EI  / 9  u  ,  1  _ 

^77* 


1  92U  _  1 

77*'* 


2.1-7 


Taking  the  limit  of  eq.  2.1-7  as  R  +  °°  and  integrating  the  resulting 
equation  twice  with  respect  to  ds  yields 


4  2 

EI  9  u  ,  9  u  -Ip 
m  as4  at2 ' m  R  ’ 


2.1-8 


where  the  integration  time  functions  are  omitted.  Equation  2.1-8  is 
the  governing  equation  for  the  flexural  response  of  a  straight  beam. 
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2.2  Solution  to  the  Homogeneous  Equation 

The  homogeneous  equation  can  be  written  as 


El  ,36w  , 

4  '  fi 
mlT  34>° 

where  w (4> , t )  represents 
of  the  form 


♦ 

3<j)  3<j> 


4 


34>23t 


_  .  l_w  =  o 
?  2  u 
3t 


2.2-1 


either  v(4>,t)  or  u(4>,t).  Assuming  a  solution 


w(4> ,t)  =  l  W  ((D)  T  (t)  2.2-2 

n=l  n  n 


separates  eq.  2.2-1  into 


. i V I  ,  .1 V  ,  / 1  -,2\  ..II  ,  ,2  •• 
Wn  +  2W  +  (1-An)  Wn  +  Xn  Wn 


=  0 


2.2-3 


and 


n 


“n  Tn 


=  0 


2.2-4 


where 


and 


.2  2  mR4 

An  =  \  E T 


4LI.  (  ) 

d4>  1  ' 


I 


dt 


(  )  ,  etc. 


2.2-5 


The  solution  to  eq.  2.2-4  is 
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2.2-6 


Since  eq.  2.2-3  is  an  ordinary  differential  equation  with  constant 
coefficients  a  solution  of  the  form 

6  r  .<f> 

-  I  Cnj  e  2 

J  * 

can  be  assumed.  Substituting  eq.  2.2-7  into  eq.  2.2-3  yields  an 
auxiliary  equation 


rjj  +  2rl  +  (1-X^)  +  xf  =  0  .  2.2-8 

n  n  n  '  n'  n 

2 

Equation  2.2-8  is  a  cubic  equation  in  rn  and  applying  the  standard 

algebraic  solution  for  cubic  equations  (Sokolnikoff  (23))  a  solution  for 
2 

rn  is  obtained.  By  considering  the  discriminant  of  the  cubic  equation, 

2 

it  was  found  that  for  0  <  Xn  <  0.113,  the  solution  to  eq.  2.2-3  can  be 
expressed  as 


W  (<|>)  =  Cn-|  cos  v-j <J)  +  Cn2  sin  v-j 4>  +  cos  V2<P 


where 
and  the  r 
for  0.113 


+  C  ^  sin  \>2$  +  Cn5  cos  v3<j>  +  Cn6  sin  v3<j)  2.2-9 

rn(l  ,2)  =  -  1V1  *  rn(3 ,4)  =  -  1v2  and  rn(5,6)  =  -  lv3 
n ( j ) ' s *  J  =  2,  ...6  are  the  roots  of  eq.  2.2-8.  Similarly 

2 

<  X  <  17.64  the  solution  is 
n 
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Wn (4>)  =  Cn-j  cos  v.j<J>  +  Cn2  sin  v^<j>  +  cos  cosh  y2<|> 


+  Cn4  cos  v24>  sinh  y24>  +  Cn5  sin  v2<t>  cosh  y2<j>  +  Cn6  sin  v2<j)  sinh  y2<j> 


where  rn(l,2)  =  -  ivl  >  rn(3,4) 
2 

For  \  >  17.64  the  solution  is 
n 


2.2-10 

±u2±iv2  and  rn(5>6)  »  ±  W2  +  iv2  . 


»„M  =  Cni  cos  v-|4>  +  Cn2  sin  v-|<|>  +  cosh  y-j<j) 


+  Cn4  sinh  y -j 4>  +  Cn5  cosh  y24>  +  Cn6  sinh  2.2-11 

where  rn(l ,2)  =  i  ivl  >  rn (3 ,4)  =  ±  h  and  rn(5,6)  *  ±  ^2  '  The 
constants  Cnj  ,  j  =  1,2,... 6  are  determined  from  the  boundary  conditions. 

Equations  2.2-9,  2.2-10  and  2.2-11  are  the  solutions  to  the  homogeneous 
eq.  2.2-3  where  Wn(4>)  represents  either  Vn(<{>)  or  Un(<j>)  depending  whether 
the  tangential  or  radial  mode  shapes  are  considered.  Before  a  particu¬ 
lar  solution  to  eq.  2.1-5  and  2.1-6  can  be  evaluated,  the  orthogonality 
condition  for  the  free  modes  of  vibrations  must  be  found. 

2.3  Orthogonality  Condition 

To  evaluate  the  orthogonality  condition,  consider  eq.  2.2-3 
where  Wn (4>)  is  replaced  by  Vn (d>) 

v!!  +  2  V^V  +  (l-xj;)  V*1  +  /  ?  =  0  . 

n  n  n  n  n  n 


2.3-1 
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Multiplying  eq.  2.3-1  through  by  Vm  c4  and  integrating  between  the 
limits  of  -a  to  a  yields 


where 


Bn  - 


a 


-a 


( v£ 1  IV_1 1 1  -2  v”v”  +  vV)  d<(>= 


n  m 


n  m  n  nr 


-  X, 


n 


a 


-a 


(vV  +  V  V  )  d* 
n  m  n  m 


Bi  =  -  vJVvI  +  +  2  V*11  Vm 

I  n  m  nm  nm  n  m 


2.3-2 


a 


- 2  v'V +  (i-x?)  vivj 

n  m  v  n7  n  m-a 


The  constant  B-j  can  be  evaluated  after  the  boundary  conditions  at 
<J)  =  +  a  are  specified.  Another  equation  similar  to  eq.  2.3-2  can  be 
obtained  by  interchanging  subscripts  m  and  n.  This  yields 


where 


B2- 


a 


-a 


(v"V!I  -  2  v’V1  +  vV)  d<J) 

m  n  m  n  m  n7 


-  x: 


m 


a 


-a 


(vV  +  VI)  d<t> 

'  m  n  m  n  r 


B?  =  CVIVn  ‘  +  Villvn1  +  2 

2  L  m  n  mn  mn  mn 


2.3-3 


a 


”  2  V*V  +  (1-A?)  vlvj 

m  n  v  n;  m  n  -a 


Subtracting  eq.  2.3-3  from  eq.  2.3-2  yields 
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B1  -  B2  ■  -  # 


a 


-a 


(vV  +  V„V J  d*. 
n  m  n  m 


2.3-4 


When  pinned-end,  clamped-end  or  free-end  boundary  conditions 
are  applied  and  B-|  and  are  evaluated,  it  is  found  that  the  left- 
hand  side  of  eq.  2.3-4  is  identically  zero.  Equation  2.3-4  can  be 
written  as 


U 


m 


o  f0,  T  T 

XT)  (VV  +  V  V  )  d<})  =  0 
n'  v  n  m  n  m 

-a 


2.3-5 


Hence,  for  m  ^  n  , 


a 


-a 


(vV  +  V „VJ  d*  ■  0  , 
n  m  n  m 


2.3-6 


and  for  m  =  n,  the  eigenvectors  are  normalized  such  that 


2a 


a 


-a 


(V¥  +  V  V)  d«* 
mm  mm 


=  1 


2.3-7 


Using  the  orthogonal  condition  as  given  by  eq.  2.3-6  and  the 
normalization  given  by  eq.  2.3-7,  the  steady  state  solution  to  eq.  2.1-5 
and  2.1-6  can  be  evaluated. 

2.4  Steady  State  Solution 

The  steady  state  solution  to  eq.  2.1-5  can  be  assumed  to  be 


v  (<!>,t)  =  l  V.(<j>)  n . (t) 

P  j=l  J  3 


2.4-1 
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where  V.  is  the  j-th  free  mode  of  vibration  and  n-  is  a  function  as  yet 

J  J 

undetermined.  Substituting  eq.  2.4-1  into  eq.  2.1-5  yields 


I  (V«  +  2  v‘V  +  V”)  n.  +  n'i  (v"-V.)] 
j=l  mir  J  J  J  J  J  J  J 


l  3Pr 

m  [d(p  V  * 


2.4-2 


Using  eq.  2.3-1,  eq.  2.4-2  can  be  reduced  to 


l  [(v"  -  VjHfij  +  4j)] 

j  =  l 


1  9PR 

j-  -  PT) 

m  v3cf>  T' 


2.4-3 


Multiplying  eq.  2.4-3  through  by  d<J>  and  integrating  between  the 
limits  of  -a  to  a  yields 


a 


I  [(fij  +  UjTlj)  J  (vjv‘  +  VjVk)  d*] 

-a 


j  =  l 


ra  Vk  8Pr 

m  d(}) 

-a 


2.4-4 


Applying  the  orthogonality  condition  as  given  by  eq.  2.3-6  uncouples 
eq.  2.4-4  and  then  using  the  normalization  as  given  by  eq.  2.3-7,  eq. 
2.4-4  can  be  written  as 


\ +  “k\ =  Dk(t) 


2.4-5 
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where 


¥*>  ■  -  £ 


a 


Vk  9PR 

m  3<t>  '  PT^  d<t' 


2.4-6 


-a 


The  solution  to  eq.  2.4-5  can  be  expressed  in  integral  form  as 


nR(t)  =  c 


k-j  sin  oo^t  +  c^2  cos  co^t 


Dk(x)  sin  u)k(t-x)  dx 


2.4-7 


where  ck-| 
arch. 


I 

and  ck2 


are  determined  from  the  initial  conditions  of  the 


The  steady  state  solution  for  the  tangential  displacements  is 
obtained  by  substituting  eq.  2.4-7  into  eq.  2.4-1.  The  steady  state 
solution  for  the  radial  displacements  can  be  obtained  in  a  manner  similar 
to  that  above,  however,  it  can  be  done  more  simply  by  using  the  in¬ 
extensibility  condition  as  given  by  eq.  1.2-12. 

Steady  state  solutions  are  presented  for  two  types  of  support 
movement.  The  motions  will  be  defined  by  "symmetrically  excited  arches" 
for  the  one  case  and  "unsymmetrically  excited  arches"  for  the  second 
case. 


2.5  Solution  for  Symmetrical  Excitation 

The  term  "symmetrical  excitation"  is  taken  to  mean  that  both 
supports  A  and  B,  in  Fig.  2.5-1,  are  vibrating  in-phase  and  in  the  plane 
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of  the  initial  curvature  of  the  arch.  Both  supports  have  displacements 
whose  direction  is  parallel  to  the  line  <|>  =  0. 


<f>=0 


FIG.  2.5-1  MODEL  FOR  SYMMETRICAL  EXCITATION 


Let  the  function  describing  the  support  movement  be  denoted  by  g(t)  where 
g(t)  is  twice  continuously  differentiable.  Under  this  type  of  excitation 
the  external  forces  P-j.  and  PR  are 


P-p  =  -  m  g(t)  sin  <j> 


2.5-1 


and 
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PR  =  -  m  g(t)  cos  <f>. 


2.5-2 


Substituting  eq.  2.5-1  and  2.5-2  into  eq.  2.4-6  and  from  eq.  2.4-7  the 
solution  for  nk(t)  is 


nk(t)  =  - 


Dk  g(x) 


sin  ook(t-x)  dx 


2.5-3 


where 


I 


a 

Vk((j))  sin  (f>  d4> 

-a 


2.5-4 


and  the  initial  conditions  can  be  disregarded  since  some  damping  is 
present  in  an  actual  arch  which  will  eliminate  any  initial  effects 
leaving  only  the  steady  state  terms.  If  g(t)  is  assumed  to  be 

g(t)  =  G  sin  fit  ,  2.5-5 


the  solution  for  rik(t)  becomes 


?  1 
G  JT  D, 


nk(t)  =  ~2 - 2  S1*n  ^ 


2.5-6 


ook  -  fi 


Substituting  eq.  2.5-6  into  eq.  2.4-1,  the  steady  state  solution  for 
the  tangential  displacement  becomes 


v  “  fi  D .  V. 

— E.  =  y  — _J — si— — 

G  A  72~27 

J=1  (Wj  -  fi  ) 


sin  fit 


2.5-7 
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and  using  eq.  1.2-12,  the  radial  displacement  is 


UJL-_  y  ^  DJ  UJ 
G  j=l  (u)j  -  fi2) 


sin  fit 


2.5-8 


where  the  subscript  p  means  steady  state  solution. 

The  terms  Vp/G  and  up/G  are  referred  to  as  magnification 
ratios  in  the  tangential  and  radial  directions  respectively. 


2.6  Solution  for  Unsymmetrical  Excitation 

With  reference  to  Fig.  2.6-1,  "unsymmetrical  excitation"  is 
taken  to  mean  that  support  A  is  stationary  and  support  B  is  vibrating 
in  the  plane  of  the  initial  arch  curvature.  The  displacement  of  B  is 
along  a  line  parallel  to  the  line  <j>  =  0  . 


FIG.  2.6-1  MODEL  FOR  UNSYMMETRICAL  EXCITATION 
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Assume  the  function  describing  the  support  movement  is  given  by  g(t) 
where  g(t)  is  twice  continuously  differentiable.  Let  £(t)  define  the 
rotation  of  the  arch  about  support  A  due  to  the  movement  of  support  B. 
Point  P  is  an  arbitrarily  chosen  point  on  the  circumference  of  the 
arch.  If  the  amplitude  of  g(t)  is  small  compared  with  the  distance  AB, 
the  rotation  of  the  arch  is  given  by 

s(t)  =  fltel 

AB 


or  by 


£(t)  = 


9(0 


2R  sin  a 


2.6-1 


since  AB  =  2R  sin  a.  The  distance  AP  can  be  expressed  as 


AP  =  2R  sin  j  U+a)  ,  2.6-2 

hence,  the  displacement  of  point  P  (represented  by  z  in  Fig.  2.6-1)  is 
given  by 


z(<J>,t)  =  AP  £(t)  .  2.6-3 

The  external  forces  Py  and  P^  are  given  by 


Py  =  -  mW  sin  ( 0+4>)  £ 


2.6-4 
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and 


where 


Substituting 


and 


Substituting 


and  from  eq. 


where 


PR  =  -  m  J?  cos  (0+4>)  'i 


2.6-5 


0  =  \  . 


for  AP,  0  and  £  in  eq.  2.6-4  and  2.6-5  yields 


pr  ■  -  rfrL  n  - cos 


2.6-6 


pr  =  -  5-frf-s  sin  (*+a)  • 


2.6-7 


eq.  2.6-6  and  2.6-7  into  eq.  2.4-6,  yields 


D  (t) _ -iii) _ 

4a  sin  a 


2.4-7  nk(t)  becomes 


a 


-a 


Vk  [1  -  2  cos  (<j>+a)]d(J) , 


nk(t)  -  ck-j  sin  wkt  ^  ck2  ^k^ 


(jO. 


g(x)  Dk  sin  wk(t-x)  dx 


2.6-8 


D,  = 


1 


k  4a  sina 


a 


-a 


V k  [1  -  2  cos  (cj>+a)]  d(J>  2.6-9 
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Substituting  the  assumed  sinusoidal  function  for  g(t)  as  given  by  eq. 
2.5-5  into  eq.  2.6-8,  ( t )  can  be  written  as 

2  ' 

G  fi  D. 

nk(t)  =  ~2 - ^  sin  fit  2.6-10 

00k  -  fi 

where  the  initial  conditions  are  disregarded  due  to  the  presence  of 
some  damping  in  a  physical  arch.  When  eq.  2.6-10  is  substituted  into 
eq.  2.4-1,  the  steady  state  solution  for  the  tangential  displacement  is 


vp(4>»t) 

"G 


n2  VjU) 

(w?  -  Q2) 

J 


sin  fit 


2.6-11 


and  using  eq.  1.2-12,  the  steady  state  solution  for  the  radial  dis¬ 
placement  is 


up(4>,t) 

"G 


00  fi^  D.  U .  ( 4> ) 

7  - r?  ^ a —  sin  fit  , 

j=l  (u?  -  fi2) 


2.6-12 


where  the  subscript  p  denotes  the  steady  state  solution.  Similarly  as 
was  defined  for  the  symmetrical  case,  the  terms  up/G  and  vp/G  are 
referred  to  as  magnification  ratios  in  the  radial  and  tangential  di¬ 
rections  respectively. 


2.7  Computer  Evaluation  of  Theoretical  Solutions 

The  theoretical  solutions  were  evaluated,  only  for  pinned-end 
boundary  conditions,  using  a  digital  computer.  The  boundary  conditions 
were  applied  to  the  homogeneous  solution  as  given  by  eq.  2.2-9,  2.2-10 
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or  by  2.2-11  and  the  constants  C  . ,  n  =  mode  number  and  j  =  1,2,... 6, 

« •  J 

were  evaluated.  The  pinned-end  boundary  conditions,  i.e.  u  =  0,  v  =  0 

and  M=  0  at  ^  =  +  a  ,  are  homogeneous  and  application  of  these  to  eq. 

2.2-9,  2.2-10  or  2.2-11  results  in  six  homogeneous  algebraic  equations. 

In  order  that  a  non-trivial  solution  for  this  set  of  equations  exist, 

the  determinant  of  the  coefficients,  which  contains  the  natural  frequency 

2 

parameter  =  wn  R  vtn/EI '  ,  must  be  zero.  Since  it  was  a  difficult 
task  to  obtain  an  exact  solution  for  An  by  expanding  and  solving  the 
six  by  six  determinant,  an  iterative  procedure  was  used.  The  procedure 
was  started  by  assuming  a  value  for  An  which  made  the  determinant 
negative  and  assuming  another  An  which  made  the  determinant  positive. 

An  intermediate  value  for  An  was  then  calculated  and  the  determinant 
evaluated  again.  This  procedure  was  continued  until  the  determinant 
was  zero  or  until  the  value  of  An  was  iterated  to  16  significant 
figures  (capacity  of  the  computer).  It  was  not  known  before  hand  to 
which  mode  the  particular  value  of  An  would  correspond;  however,  the 
mode  was  identified  when  the  mode  shape  was  evaluated.  In  the  iterative 
process,  it  was  possible  to  have  the  resulting  determinant  approach 
zero  for  the  lower  modes,  however  for  the  higher  modes  a  small  change 
in  the  value  of  An  (i.e.  a  change  in  the  16th  significant  figure) 
produced  large  changes  in  the  value  of  the  determinant.  For  example, 
for  an  arch  having  a  *  100°,  a  change  by  one  in  the  16th  significant 
figure  for  Ag  caused  the  value  of  the  determinant  to  change  from  *  10^° 
to  ~  -  1010.  However,  for  the  same  arch  it  was  possible  to  iterate  for 
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_1  o 

A-|  such  that  the  value  of  the  determinant  was  ~  10  .  The  value  of 

the  determinant  was  not  zero  for  the  higher  modes,  but  because  of  the 
large  oscillations  in  the  value  of  the  determinant  for  small  changes 
in  An,  it  was  decided  that  the  iterative  method  evaluated  the  An's  with 
sufficient  accuracy. 

To  aid  in  selecting  the  initial  values  of  An  for  the  iterative 
process  for  arches  with  various  a  values,  a  plot.  Fig.  2.7-1,  was  made 
of  An  versus  a  for  the  first  nine  modes  of  free  vibrations  for  pinned- 
end  boundary  conditions.  A  similar  plot  for  fixed-end  boundary  condi¬ 
tions  is  shown  in  Fig.  2.7-2. 

A  comparison  can  be  made  of  the  values  of  An  obtained  by  the 
present  theory  and  those  obtained  by  Den  Hartog  (2),  Eppink  and 
Veletsos  (10)  and  Nelson  (17,21)  who  used  approximate  methods.  As 
mentioned  previously,  an  error  was  noted  in  the  formula  Den  Hartog 
derived  for  A-|  when  displacements  in  the  form  of  a  polynomial  were  as¬ 
sumed.  In  Table  2.7-1,  results  are  presented  for  A-j  as  obtained  using 
the  present  theory,  using  Den  Hartog' s  tabulated  values,  using  the 
sinusoidal  displacements  and  the  equation  derived  by  Den  Hartog,  and 
using  the  polynomial  displacements  assumed  by  Den  Hartog  and  the  cor¬ 
rected  formula  (i.e.  A-j  =  l/a^[(a^  +  79.20  c?  +  1584)/0  +  0.07769  a^)]^ 
for  A-j.  As  expected,  the  values  Den  Hartog  obtained  using  the  Rayleigh 
quotient  are  higher  or  equal  to  the  values  given  by  the  present  theory. 

Eppink  and  Veletsos  presented  values  only  for  one  arch  which 
had  a  =  43.6  degrees.  Results  were  presented  when  the  arch  was  approxi¬ 
mated  by  different  number  of  bars.  On  the  basis  of  this  data,  arches 
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FIG.  2.7-1. DIMENSIONLESS  NATURAL  FREQUENCY; 
PINNED-END  BOUNDARY  CONDITION 
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FIG.  2.7-2  DIMENSIONLESS  NATURAL  FREQUENCY; 

FIXED-END  BOUNDARY  CONDITION 
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TABLE  2.7-1  COMPARISON  WITH  DEN  HARTOG'S  RESULTS 
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approximated  by  less  bars  gave  lower  values  for  the  natural  frequencies 
than  those  approximated  by  a  larger  number  of  bars,  which  indicates 
that  this  method  yields  lower  bounds  for  the  natural  frequency.  On  the 
basis  of  a  20  bar  approximation,  the  values  for  A  ,  n  =  1,2,... 9  and 
the  values  obtained  from  the  present  theory  are  shown  in  Table  2.7-2 
for  an  arch  with  a  =  43.6  degrees. 


xi 

a2 

A3 

A4 

A5 

*6 

A7 

A8 

X9 

*Eppink  & 
Veletsos’s 

14.8 

36.4 

65.8 

103.1 

150.9 

205.1 

269.6 

340.1 

420.6 

Results 

Present 

Theory 

14.8 

34.6 

65.8 

102.9 

151.0 

205.1 

270.4 

341.5 

423.8 

*  Results  for  a  20-bar  approximation 


TABLE  2.7-2  COMPARISON  WITH  EPPINK  AND  VELETSOS'S  RESULTS 

From  Table  2.7-2,  Eppink  and  Veletsos  obtained  the  same  or 
lower  values  for  An  in  all  cases  except  when  n  =  2  and  n  =  4.  No  ex¬ 
planation,  other  than  round  off  error,  can  be  given  for  this. 

A  comparison  with  the  values  Nelson  calculated  for  the  first 
four  modes  is  shown  in  Table  2.7-3.  It  should  be  noted  that  Nelson's 
values  were  obtained  from  a  graph  and  may  not  be  accurate  to  the  number 
of  significant  figures  listed.  As  expected,  Nelson  obtained  larger 
values  for  the  An's  using  the  Rayleigh-Ritz  technique  in  conjunction 


with  Lagrangian  multipliers  than  those  obtained  from  the  closed  form 
solution. 
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\a,Deg 

X 

* 

0 

45 

90 

135 

180 

N 

00 

13.8 

2.27 

0.474 

0 

A1 

T 

00 

13.764 

2.267 

0.474 

0 

x2 

N 

00 

32.5 

6.93 

2.39 

0.901 

T 

00 

32.404 

6.923 

2.366 

0.901 

N 

00 

62.0 

14.1 

5.44 

2.46 

A3 

T 

00 

61.673 

13.978 

5.349 

2.447 

N 

00 

97.8 

23.3 

9.32 

4.61 

A4 

T 

00 

96.446 

22.820 

9.267 

4.597 

*N  -  Nelson's  values,  obtained  from  graph 
T  -  Values  obtained  from  present  theory 


TABLE  2.7-3  COMPARISON  WITH  NELSON'S  RESULTS 


Since  the  agreement  between  the  values  of  XR  for  the  present 
theory  and  those  of  previous  investigators  is  favourable,  it  can  be  as¬ 
sumed  that  the  iterative  method  used  in  evaluating  An  is  sufficiently 
accurate.  After  the  values  of  An  were  evaluated,  Cnj.'s  were  calculated 
using  the  "least  squares  method".  For  example,  suppose  the  original 
system  of  equations  were  defined  by 


J 


Au 


■ 


0 


2.7-1 
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where  n  equals  the  mode  number  and  the  A.  .'s  are  constants  which  are 

•  0 

evaluated  when  An's  are  known.  Dividing  eq.  2.7-1  by  CmL,  4  0,  1  £  k  £  6 


nk 


results  in  five  unknown  ratios  and  6  nonhomogeneous  equations  which 


are  inconsistent  since  the  exact  value  of  A„  is  not  known.  The  re- 

n 

suiting  equation  is 


6  C  . 

I  [f^-  A..  (A  )  +  (A  )]  =  0  .  2.7-2 

j,i  =  l  Unk  13  n  1K  n 
j^k 


Equation  2.7-2  was  reduced  to  five  equations  by  defining  a  new  function 


and  then  minimizing  e  with  respect  to  each  of  the  five  ratios  Cn j/Cnk . 

The  ratios  were  then  evaluated  by  solving  simultaneously  the  resulting 
five  equations.  Equation  2.3-7  was  used  to  solve  for  C 

When  the  An's,  Un's  and  Vp's  were  known,  the  steady  state 
solution  for  radial  displacements,  as  given  by  eq.  2.5-8  for  symmetrical 
excitation,  and  by  eq.  2.6-12  for  unsymmetrical  excitation  was  evalu¬ 
ated.  For  pinned-end  boundary  conditions,  it  was  observed  that  if 
steady  state  solutions  were  required  near  the  fourth  mode  it  was  suffi¬ 
cient  to  use  eight  terms  in  the  series  to  obtain  an  accuracy  to  three 
significant  figures.  Better  accuracy  was  obtained  for  forcing  frequencies 
below  the  fourth  natural  frequency  when  an  eight  term  series  was  used. 


2.8  Closure  to  Chapter  II 
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A  theory  has  been  presented  which  can  be  used  to  evaluate  the 
dynamic  response  of  thin  circular  arches  subjected  to  inertial  loads 
caused  by  cyclic  support  movement.  The  arch  was  assumed  to  be  inex- 
tensible,  radial  displacements  were  small  such  that  R  +  u  ~  R,  effects 
of  damping  on  the  steady  state  solution  were  assumed  to  be  negligible, 
displacements  were  due  to  flexural  effects  only  and  the  displacements 
were  in  the  initial  plane  of  the  arch  curvature.  The  orthogonality 
condition  for  the  free  modes  of  vibration  was  found  which  was  used  in 
the  evaluation  of  the  steady  state  solution.  Steady  state  solutions 
were  evaluated  for  both  symmetrical  and  unsymmetrical  excitations.  For 
forcing  frequencies  near  the  fourth  mode,  an  eight- term  series  in  terms 
of  the  free  modes  of  vibration,  yielded  the  steady  state  solution 
accurate  to  three  significant  figures. 

Since  experimental  results  will  only  be  presented  for  pinned- 
end  arches,  the  theory  has  been  evaluated  in  detail  for  these  boundary 
conditions.  Details  of  the  experimental  procedure  are  given  in  Chapter 


III. 
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Chapter  III 

Experimental  Procedure 

3. 1  Description  of  Testing  Equipment  and  Arches 

The  arches  were  excited  by  an  Unhol tz-Dickie,  Model  No.  100 
Electrodynamic  Shaker  Table,  capable  of  generating  a  maximum  force  of 
225  pounds.  The  magnitude  of  the  sinusoidal  displacement  and  the  fre¬ 
quency  of  excitation  were  controlled  by  a  Bruel  &  Kjaer  Model  No.  11019 
Automatic  Vibration  Exciter  Control.  With  this  equipment  any  desired 
frequency  and  displacement  could  be  obtained,  within  the  limitations 
of  the  apparatus.  The  frequency  was  measured  by  a  General  Radio  Model 
No.  1150-AP  Digital  Frequency  Meter  having  an  accuracy  of  plus  or  minus 
one  count.  During  testing,  the  frequency  was  recorded  continuously  on 
a  General  Radio  Data  Printer  Model  No.  1137-A  in  order  to  determine 
whether  drifting  of  the  frequency  occurred  during  a  test.  The  amplitude 
of  displacement  for  the  support  and  that  of  the  arches  were  measured 
by  micrometer  probes,  shown  in  Fig.  3.1-1.  The  probes  were  fitted  with 
spring  loaded  tips  to  reduce  any  disturbances  that  might  be  induced 
when  contact  was  made  with  the  arch  or  support.  The  probe  measuring 
the  radial  arch  displacement  was  mounted  on  a  movable  arm  permitting 
the  rotation  of  the  probe  into  any  desired  position.  A  simple  electrical 
circuit  was  used  to  indicate  when  the  probe  made  contact  with  the  arch. 

A  photographic  view  of  the  testing  equipment  is  shown  in  Fig.  3.1-2. 
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FIG.  3.1-1  MICROMETER  PROBES 


FIG.  3.1-2  GENERAL  VIEW  OF  TESTING  EQUIPMENT 
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The  arches  were  rolled  from  strips  cut  from  6061-T6  aluminum 
sheets  having  a  modulus  of  elasticity  of  10^  psi.  The  arches  were  bolted 
to  one-half  inch  diameter  steel  shafts  with  bearing  supports  and  in 
order  to  reduce  the  weight  of  the  model,  were  mounted  on  a  3"  x  1-1/2"  x 
1/4"  aluminum  channel.  To  reduce  the  effects  of  support  flexure,  the 
length  of  the  supporting  channel  was  chosen  to  be  11-inches. 

To  locate  the  points  of  zero  displacement  on  the  arch,  a 
Chadwick-Helmuth  Model  No.  126A  Stroboscope  synchronized  with  the  support 
displacement  and  equipped  with  a  Chadwick-Helmuth  Slip-Sync  Model  No. 

105  AR  Automatic  Phase  Shift  was  used. 

3.2  Testing  Procedure 

Any  arbitrary  point  on  the  arch  has  both  radial  and  tangential 
displacements  which  are  related  by  eq.  1.2-2.  Because  of  this  relation¬ 
ship,  experimental  evaluation  of  only  one  displacement  was  considered. 
Since  radial  displacements  were  easier  to  measure  the  theory  is  evalu¬ 
ated  on  this  basis. 

It  was  decided  to  use  pinned-end  supports  for  the  experimental 
tests  because  previous  experience  with  vibrating  cantilever  beams 
indicated  that  fixed-end  conditions  were  difficult  to  attain  [see  Rae 
(24)]. 

To  determine  the  accuracy  of  the  micrometer  measuring  probe, 
a  test  was  conducted  comparing  displacements  as  measured  by  the  probe 
against  those  measured  by  a  Gaertner  Cathetometer  which  is  an  optical 
measuring  device.  The  comparison  is  shown  in  Fig.  3.2-1.  In  view  of 
the  small  difference,  the  micrometer  probe  was  chosen  as  a  reliable 
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CATHETOMETER  REAOING  (INCHES) 

FIG.  3.2-1  COMPARISON  OF  DISPLACEMENT  READINGS 
BETNEEN  MICROMETER  PROBE  AND  CATHETOMETER 
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means  for  measuring  deflections.  The  precision  of  the  probes  was  found 
to  be  approximately  +0.005  inches. 

In  order  to  determine  what  arch  parameters  would  yield  the 
most  reliable  experimental  data,  a  series  of  arches  with  various  R/H 
and  a  values  was  tested.  The  arch  behavior  was  observed  under  both 
symmetrical  and  unsymmetrical  excitation.  It  was  concluded  that  arches 
with  a  *  140  degrees  and  43  £  R/H  £  128  were  very  susceptible  to  dis¬ 
turbances  induced  by  the  measuring  probe.  When  arches  with  a  <  90 
degrees  and  R/H  *  130  were  tested,  they  were  very  stiff  and  it  was 
impossible  to  obtain  a  measurable  displacement  response  at  frequencies 
other  than  the  first  resonant  frequency.  Because  of  the  force  limitation 
of  the  shaker  table  and  because  of  the  resolution  of  the  displacement 
measurement,  it  was  found  that  the  best  operating  range  of  the  shaker 
table  for  these  tests  was  between  5  and  125  H  .  On  this  basis,  it  was 
decided  to  build  arches  which  had  the  first  resonant  frequency  above 
5Hz  and  at  least  the  second  below  125  H  .  With  these  frequency  limitations, 
a  forced  response  could  then  be  obtained  for  frequencies  up  to  the 
second  mode,  although  for  some  arches  a  forced  response  was  obtained 
for  frequencies  up  to  the  fourth  mode.  To  meet  the  above  limitations, 
it  was  decided  the  best  range  of  arches  to  test  was  90°  £  a  £  125°  which 
required  the  radius  to  be  5.01  £  R  £  6.00  since  the  distance  between 
supports  was  fixed  at  10.02  inches  for  all  tests.  The  condition  of 
inextensibility  i.e.  (R/k)^  »  \ ^  was  satisfied  if  thicknesses  of  0.032 
and  0.042  inches  and  a  width  of  one  inch  were  used. 

It  was  also  observed  that  arches  having  large  variations  in 
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radius  had  experimental  values  considerably  higher  than  those  predicted 
by  theory.  Because  of  the  higher  experimental  values,  this  would 
indicate  that  some  unknown  parameter  is  highly  sensitive  at  this  fre¬ 
quency  and  is  exciting  the  arch  amplitude  beyond  the  theoretical  value. 
The  best  circular  arches  were  obtained  by  using  a  Di-Arco,  No.  2A 
precision  three-roll  roller.  The  arches  were  checked  by  tracing  the 
arch  outline  on  paper  and  then  drawing  in  the  best-fit  radius.  Measure¬ 
ments  were  taken  and  a  radius,  R  +  ^-H,  was  assigned.  The  theoretical 
solutions  were  based  on  the  radius  (R  +  Jr  H)  which  most  closely  matched 
the  theoretical  and  experimental  values  of  the  resonant  frequency  for 
the  first  four  modes.  Any  arch  having  a  radius  variation  greater  than 
one-half  the  thickness  was  rejected.  The  values  of  a  for  the  arches 
were  calculated  on  the  basis  of  the  best-fit  radius  and  by  using  the 
distance  between  the  arch  supports  of  10.02  inches. 

As  noted  in  Chapter  II,  under  the  assumptions  used  in  the 
development  of  the  theory,  some  damping  is  present  in  the  system,  how¬ 
ever,  it  was  assumed  to  be  small  and  was  neglected  in  the  theoretical 
analysis.  To  reduce  the  effects  of  damping  in  the  experimental  tests, 
arches  were  excited  at  frequencies  sufficiently  far  from  resonant 
conditions  since  the  effect  of  damping  is  most  pronounced  at  resonant 
conditions.  To  meet  this  condition,  it  was  decided  to  excite  arches 
at  frequencies  such  that  the  maximum  magnification  ratio  ( U/G)  was 
kept  below  the  value  of  four. 

To  check  the  effects  of  damping  in  the  system,  one  LVDT 
(linear  variable  differential  transformer)  probe  was  attached  to  the 
crown  of  the  arch  and  another  LVDT  probe  was  attached  to  the  support. 
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Under  symmetrical  excitation  (to  ensure  no  tangential  displacement  at 
the  crown)  a  simultaneous  trace  of  the  output  from  the  two  LVDT's  was 
recorded,  using  a  Brush  Oscillograph  Model  No.  13-5335-00,  at  various 
frequencies.  As  shown  in  Fig.  3.2-2,  the  two  outputs  are  nearly  in- 
phase  or  180  degrees  out-of-phase.  On  the  basis  of  these  phase  angles, 
it  was  concluded  that  damping  was  small  and  could  be  omitted  from  the 
theoretical  analysis.  Although  a  frequency  trace  was  obtained  only  for 
one  point  on  the  arch,  other  points  were  assumed  to  behave  similarly. 

By  comparing  displacements  obtained  from  taking  readings  on 
the  outside  and  inside  circumference  of  the  arch,  a  check  was  made  on 
the  zero  shift  of  the  static  position  of  the  arch.  Since  the  differences 
in  the  readings  were  in  the  range  of  the  precision  of  the  micrometer 
probe,  it  was  decided  that  sufficiently  reliable  results  would  be  ob¬ 
tained  by  taking  readings  on  the  outside  circumference  only. 

Since  the  arch  vibrated  about  its  static  position,  the 
amplitudes  of  the  radial  displacements  were  obtained  by  taking  the 
difference  in  the  readings  between  the  dynamic  and  static  reading. 

For  example,  consider  an  arbitrary  point  P,  shown  by  the  static  position 
A  in  Fig.  3.2-3,  for  which  a  reading  was  taken.  Point  P  is  moving  con¬ 
tinuously  and  reaches  its  maximum  displacement  when  P  is  at  P'  in 
position  B,  for  which  a  dynamic  reading  was  taken.  The  difference  in 
the  two  readings  is  U^»  where  represents  the  absolute  radial  dis¬ 
placement.  In  order  to  eliminate  the  error  e  which  would  occur  if  both 
the  static  and  dynamic  readings  were  taken  with  the  probe  at  the  same 
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ARCH  (<f>=0) 

SUPPORT 


<  oo2  (fl  =  1  9.3  H  ) 


ARCH  (<J>=0) 


SUPPORT 


co2  <  <  034  [o,  =  60,4  H  ) 


ARCH  (4>=0) 


SUPPORT 


fi  >  co4  (fi  =  160.0  H  ) 


FIG.  3.2-2  DISPLACEMENT  TRACES  FOR  ARCH  AND  SUPPORT 
R/H  =  133.3  and  a  =  116.5  DEGREES 
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angular  position,  it  was  necessary  to  move  the  probe  to  position  P'  when 
taking  the  dynamic  reading.  Generally,  the  tangential  displacement  for 
P  was  less  than  0.15  inches,  giving  a  maximum  angular  difference  in 
the  two  positions  of  less  than  two  degrees  for  an  arch  with  a  five- 
inch  radius.  Experimental  data  was  taken  at  two-inch  intervals  along 
the  outside  circumference. 

For  each  frequency  a  series  of  four  readings,  each  with  a 
different  support  displacement,  was  taken  for  each  position  on  the 
arch.  The  support  displacements  were  divided  into  four  equally  spaced 
increments  between  the  two  extremes  used  for  the  particular  test.  An 
average  value  of  absolute  arch  displacement  divided  by  the  support  dis¬ 
placement  (i.e.  U^/G)  was  calculated  from  the  four  readings  for  each 
posi tion. 

Probe  position  for 


FIG.  3.2-3  POSITION  OF  MICROMETER  PROBE  DURING  TEST 
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When  arches  were  excited  near  a  resonant  frequency,  it  was 
observed  that  the  amplitude  of  the  arch  displacement  was  highly  sensitive 
to  slight  variations  in  the  forcing  frequency.  (This  is  another  indi¬ 
cation  that  damping  was  small.)  If  UA/G  <4,  the  forcing  frequency 
would  be  sufficiently  removed  from  the  resonant  frequency  yet  be  close 
enough  to  the  resonant  mode  ensuring  a  measurable  arch  amplitude. 
Generally,  the  maximum  value  of  U^/G  was  kept  between  one  and  two. 

It  was  also  observed  for  some  arches  excited  beyond  the  fourth  resonant 
frequency,  out-of-plane  modes  were  excited.  These  out-of-plane  vi¬ 
brations  occurred  at  frequencies  approximately  sixteen  times  the  fre¬ 
quency  of  the  first  flexural  mode  which  agrees  with  the  observations 
made  by  Volterra  and  Morell  (15).  These  out-of-plane  vibrations  are 
excited  because  of  the  difficulty  in  aligning  the  plane  of  displace¬ 
ment  with  the  plane  of  the  arch. 

3.3  Presentation  of  Theoretical  and  Experimental  Results 

Results  are  presented  for  four  arches  with  R/H  values  of 
121.2,  133.3,  157.8  and  179.4  and  a  values  of  100.2,  116.5,  97.2  and 
119.2  degrees  respectively.  For  each  arch,  the  first  four  mode  shapes 
are  presented,  followed  by  the  results  for  symmetrical  excitation,  and 
then  by  the  results  for  unsymmetri cal  excitation.  For  both  symmetrical 
and  unsymmetri cal  excitations,  results  are  presented  for  four  different 
forcing  frequencies. 

Since  damping  was  small,  the  experimental  resonant  frequencies 
were  determined  by  adjusting  the  forcing  frequency  until  maximum  de- 
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flections  occurred  in  the  arch.  This  frequency  could  be  obtained 
within  +  0.3  H .  To  obtain  the  mode  shapes  of  free  vibration,  the 
arch  amplitudes  were  measured  as  the  arch  vibrated  at  the  resonant 
frequency.  Figure  3.3-1  shows  a  photographic  view  of  the  first  four 
modes  for  various  positions  of  the  arch  as  it  moved  through  its  cycle. 

At  resonant  conditions,  a  support  displacement  of  less  than  0.005 
inches  was  sufficient  to  have  arch  amplitudes  in  the  order  of  0.2  inches 
giving  a  magnification  factor  with  values  greater  than  40.  Since  it 
was  possible  to  obtain  large  magnification  ratios  at  resonant  frequencies, 
this  indicated  that  at  a  resonant  frequency,  the  arch  was  essentially 
under  the  condition  of  free  vibration.  Hence,  under  these  conditions, 
sufficiently  accurate  mode  shapes  could  be  measured.  For  the  unsymmetri cal 
modes,  the  experimental  data  was  normalized  such  that  the  largest  ex¬ 
perimental  displacement  coincided  with  the  displacement  on  the  theoretical 
curve  for  the  same  point.  For  the  symmetrical  modes,  the  experimental 
data  was  normalized  with  respect  to  the  point  at  <j>  =  0. 

The  theoretical  solution  as  given  by  eq.  2.5-8  for  symmetrical 
excitation  and  by  eq.  2.6-12  for  unsymmetri cal  excitation  yields  the 
radial  displacement  due  to  arch  deformation  only,  i.e.  the  effect  of 
support  movement  on  the  arch  displacement  is  not  considered.  The  experi¬ 
mental  values  are  absolute  displacements  which  requires  a  modification 
of  the  theoretical  solution  to  account  for  the  support  movement  before 
a  comparison  between  theory  and  experiment  can  be  made.  The  modified 
theoretical  solution  is 
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FIG.  3.3-1  FIRST  FOUR  FREE  MODES;  R/H  =  163.4  AND  a  =  106.7  DEGREES 


THEORY  23.7  HZ  -  THEORY  162.1  HZ 

O  EXP.  24.1  HZ  O  EXP.  157.0  HZ 
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FIG.  3.3-2  NORMALIZED  RADIAL  DISPLACEMENT,  R/H=121 
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FIG.  3.3-3  MRGNIFICRTION  RATIO;  SYMMETRICALLY  FORCED,  R/H=121. 
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FIG.  3.3-4  MAGNIFICATION  RATIO;  UNSYMMETRICALLY  FORCED,  R/H=121. 
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FIG.  3.3-5  NORMALIZED  RADIAL  DISPLACEMENT,  R/H=133.3 
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FIG.  3.3-6  MflGNIFICRTION  RATIO;  SYMMETRICALLY  FORCED,  R/FN133.3 
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FIG.  3.3-7  MAGNIFICATION  RATIO;  UNSYMMETRICALLY  FORCED,  R/H=133.3 
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FIG.  3.3-8  NORMALIZED  RADIAL  DISPLACEMENT,  R/H=157.8 
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FIG.  3.3-9  MAGNIFICATION  RATIO;  SYMMETRICALLY  FORCED.  R/H=157.8 
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FIG.  3.3-10  MAGNIFICATION  RATIO;  UNSYMMETRICALLY  FORCED,  R/H=157.8 
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FIG.  3.3-11  NORMALIZED  RADIAL  DISPLACEMENT,  R/H=179.4 
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FIG.  3.3-12  MflGNIFICflTION  RATIO;  SYMMETRICALLY  FORCED,  R/H=179.4 
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FIG.  3.3-13  MflGNIFICflTION  RATIO;  UNSYMMETRICALLY  FORCED,  R/H=179.4 
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for  symmetrical  excitation,  and 


ft  DjUj(cj)) 

2  2 

oj  .  -  fi 

J 


sin((j)+a) 

2sin  a 


■]  sin  fit 


3.3-2 


for  unsymmetrical  excitation,  where  denotes  absolute  radial  displace¬ 
ment  and  G  is  the  amplitude  of  support  displacement. 

A  comparison  between  theory  and  experiment  for  the  mode  shapes 
and  forced  excitations  is  shown  in  Fig.  3.3-2  through  to  3.3-13. 


3.4  Discussion  of  Theoretical  and  Experimental  Results 

A  summary  is  given  in  Table  3.4-1  comparing  the  theoretical 
and  experimental  values  for  the  resonant  frequencies  together  with  the 
percentage  difference.  It  can  be  seen  that  good  agreement  was  obtained. 
The  experimental  values  of  the  resonant  frequencies  are  generally  higher 
than  the  theoretical  ones  for  modes  one  and  two  and  lower  for  modes 
three  and  four.  A  possible  cause  for  this  is  friction  acting  at  the 
bearing  supports  forcing  the  experimental  values  to  be  higher  at  the 
lower  modes.  However,  at  the  higher  modes,  it  is  possible  that  the  in¬ 
ertia  of  the  support  shafts  overcomes  the  frictional  effects  resulting 
in  lower  experimental  values  than  those  predicted  by  theory.  Good 
agreement  between  theory  and  experiment  for  all  mode  shapes  was  ob¬ 
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162.1 
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B.  =  121  2 
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78.4 
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+  1.7 

-  0.3 

-  3.1 
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11.3 

43.8 

93.8 

158.3 

jj-  =  133.3 
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11.8 

44.3 

94.0 

154.5 
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+  4.4 

+  1.1 
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-  2.4 
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00 
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21 .2 

67.4 

131.5 

206.7 

D 

+  5.5 

+  2.6 

-  2.1 

-  6.6 

T 

7.4 

29.9 

64.4 

109.1 

jj-  =  179.4 

E 

7.8 

30.4 

64.4 

106.3 

D 

+  5.4 

+  1.7 

0 

-  2.6 

* 

T  Theoretical  Frequency,  Hz 
E  Experimental  Frequency,  Hz 
D  Percent  Difference 


TABLE  3.4-1  COMPARISON  OF  THEORETICAL  AND 
EXPERIMENTAL  RESONANT  FREQUENCIES 
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For  purposes  of  discussion,  a  position  with  the  largest  dis¬ 
placement  was  chosen  for  each  arch  at  each  frequency  tested  and  a 
summary  of  those  points  is  shown  in  Table  3.4-2  for  symmetrical  excitation 
and  in  Table  3.4-3  for  unsymmetrical  excitation.  It  is  not  intehded  to 
mean  that  this  position  is  representative  of  the  remaining  positions 
on  the  arch,  but  rather  as  a  position  having  the  most  precise  experi¬ 
mental  value. 

For  forced  frequencies,  good  agreement  between  theory  and  experi¬ 
ment  was  obtained  for  all  tests  except  for  four  frequencies  as  shown 
in  Fig.  3.3-6  at  53.0  and  80.2  Hz>  Fig.  3.3-7  at  51.9  Hz  and  Fig.  3.3-12 
at  43.0  Hz<  It  is  possible  to  account  for  a  maximum  error  of  20  per¬ 
cent  by  considering  the  precision  of  the  micrometer  probes,  errors  in 
measuring  the  arch  radii  and  errors  due  to  variations  in  forcing  fre¬ 
quencies.  Since  at  the  forcing  frequencies  indicated  above,  the  dif¬ 
ference  in  amplitude  between  theory  and  experiment  was  greater  than 
20  percent,  this  error  is  not  due  to  experimental  error  alone.  As 
noted  previously,  arch  displacements  were  very  sensitive  to  variations 
in  radius  at  certain  frequencies  and  these  sensitivities  were  more  pro¬ 
nounced  for  arches  with  larger  a  values.  This  agrees  with  the  results 
shown  in  Table  3.4-2  and  3.4-3.  The  particular  arches  tested  were  not 
perfectly  circular,  hence  it  is  possible  that  the  four  frequencies 
yielding  poor  agreement  were  in  the  frequency  range  where  the  arch 
sensitivities  were  the  greatest.  Agreement,  within  experimental  error, 
was  obtained  for  the  remaining  28  forcing  frequencies  presented. 


* 
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Forcing 
^\^Freq . 
Arch 

Parameters 

* 

ft-j 

ft^ 

ft3 

^4 

1 

20.1 

58.1 

'  85.7 

100.2 

R/H  =  121 .2 

2 

0 

0 

0 

-  67.5 

a  =  100.2° 

3 

0.140-0.066 

0.085-0.040 

0.049-0.020 

0.055-0.026 

R  =  5.09" 

4 

0.144-0.067 

0.160-0.075 

0.156-0.072 

0.087-0.044 

H  =  0.042 

5 

1  .0% 

0.8% 

6.3% 

3.5% 

6 

<5% 

<5% 

12% 

<5% 

1 

27.7 

**53.0 

**80.2 

143.3 

R/H  =  133.3 

2 

0 

81  .9 

-  61  .4 

-  40.9 

a  =  116.5° 

3 

0.130-0.040 

0.084-0.040 

0.086-0.040 

0.025-0.011 

R  =  5.60" 

4 

0.201-0.064 

0.163-0.083 

0.101-0.047 

0.068-0.031 

H  =  0.042" 

5 

2.4% 

3.2% 

0.5% 

3.8% 

6 

10% 

42% 

25% 

10% 

1 

19.8 

35.4 

52.1 

85.0 

R/H  =  157.8 

2 

0 

0 

0 

-  68.1 

a  =  97.2° 

3 

0.138-0.050 

0.099-0.030 

0.106-0.040 

0.079-0.035 

R  =  5.05" 

4 

0.145-0.054 

0.130-0.039 

0.217-0.089 

0.120-0.059 

H  =  0.032" 

5 

1.6% 

1  .5% 

4.7% 

5.4% 

6 

<5% 

<5% 

<5% 

<5% 

1 

22.9 

**43.0 

80.0 

122.5 

R/H  =  179.4 

2 

0 

-  79.9 

59.9 

0 

a  =  119.2° 

3 

0.097-0.025 

0.112-0.026 

0.083-0.020 

0.038-0.014 

R  =  5.74" 

4 

0.212-0.055 

0.145-0.037 

0.107-0.026 

0.070-0.027 

H  =  0.032 

5 

1  .0% 

5.4% 

0.8% 

5.1% 

6 

20% 

60% 

8% 

17% 

*  1  ft,  Forcing  Frequency,  H2 

2  4> ,  Angular  Position  of  Point  Considered,  Degrees 

3  G,  Range  of  Support  Displacement  used  in  Test,  Inches 

4  U-,  Range  of  Measured  Absolute  Radial  Displacements  for  the  Point, 

A  Inches 

5  Maximum  Variation  of  UA/G  From  Average 

6  Difference  Between  Theory  and  Experiment  for  the  Point 
**  Large  Difference  Between  Theory  and  Experiment 


TABLE  3.4-2  SUMMARY  OF  RESULTS  FOR  SYMMETRICAL  EXCITATION 
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Forcing 
^\Freq . 
Arch 

Parameters 

* 

ft-j 

^2 

ft3 

^4 

1 

18.7 

34.9 

66.7 

92.0 

R/H  =  121  .2 

2 

-  45.0 

45.0 

0 

67.5 

a  =  100.2° 

3 

0.123-0.050 

0.152-0.030 

0.101-0.040 

0.080-0.034 

R  =  5.09" 

4 

0.115-0.043 

0.148-0.034 

0.144-0.055 

0.095-0.047 

H  =  0.042" 

5 

4.8% 

11.1% 

3.0% 

7.7% 

6 

10% 

<5% 

<5% 

<5% 

1 

30.1 

**51 .9 

84.7 

109.3 

R/H  =  133.3 

2 

20.5 

61.4 

40.9 

81.9 

a  =  116.5° 

3 

0.111-0.050 

0.105-0.041 

0.084-0.040 

0.055-0.029 

R  =  5.60" 

4 

0.109-0.049 

0.128-0.055 

0.135-0.063 

0.063-0.038 

H  =  0.042 

5 

0.6% 

5.1% 

CM 

• 

8.4% 

6 

<5% 

18% 

9% 

<5% 

1 

28.0 

54.6 

77.3 

123.2 

R/H  =  157.8 

2 

45.4 

0 

68.1 

22.7 

a  =  97.2° 

3 

0.135-0.060 

0.114-0.045 

0.098-0.040 

0.043-0.025 

R  =  5.05" 

4 

0.153-0.073 

0.132-0.052 

0.114-0.054 

0.064-0.035 

H  =  0.032 

5 

4.0% 

1.1% 

8.2% 

2.8% 

6 

10% 

12% 

<5% 

<5% 

1 

20.0 

39.8 

70.1 

91  .9 

R/H  =  179.4 

2 

20.0 

59.9 

-  39.9 

-  39.9 

a  =  119.2° 

3 

0.122-0.050 

0.112-0.050 

0.082-0.035 

0.079  -  0.035 

R  =  5.74" 

4 

0.116-0.050 

0.107-0.051 

0.144-0.070 

0.093-0.039 

H  =  0.032 

5 

2.9% 

3.7% 

5.7% 

5.4% 

6 

<5% 

8% 

<5% 

15% 

■k 

1  ft,  Forcing  Frequency,  Hz 

2  <|>,  Angular  Position  of  Point  Considered,  Degrees 

3  G,  Range  of  Support  Displacement  Used  in  Test,  Inches 

4  U.,  Range  of  Measured  Absolute  Radial  Displacement  for  the  Point 

A  Inches 

5  Maximum  Variation  of  U^/G  From  Average 

6  Difference  Between  Theory  and  Experiment  for  the  Point 
**  Large  Difference  Between  Theory  and  Experiment 


TABLE  3.4-3  SUMMARY  OF  RESULTS  FOR  UNSYMMETRICAL  EXCITATION 
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Other  arches  in  the  range  121.2  <  Jj-  <  179.4  and  97  .2°<  a  < 
119.2°were  also  tested,  however  the  results  have  not  been  presented  as  the 
agreement  between  theory  and  experiment  was  equally  as  good  as  those 
presented  here.  It  was  believed  that  the  four  arches  chosen  for  dis¬ 
cussion  adequately  indicated  the  effects  of  changing  arch  parameters 
( i .e.  R/H  and  a) . 

3.5  Closure  to  Chapter  III 

Experimental  verification  has  been  provided  for  the  theory 
predicting  the  dynamic  response  of  circular  thin  elastic  arches  under 
inplane  flexural  deformations.  Although  only  arches  with  simply  sup¬ 
ported  boundary  conditions  having  121  .2  £  ~  _<  179.4,  97  ,2°<  a  £  1 19.2° 
and  thicknesses  of  0.032  and  0.042  inches  were  tested,  the  theory  should 
be  valid  for  arches  having  other  arch  parameters  and  boundary  conditions. 
Considerable  care  would  be  required  In  making  the  test  arches  for  a  > 

125°  if  reliable  results  are  to  be  obtained.  To  test  arches  with  a  < 

90°  a  shaker  table  with  a  force  output  of  greater  than  225  lbs.  would 
be  required. 

Generally,  better  agreement  between  theory  and  experiment  was 
obtained  for  unsymmetrical  than  for  symmetrical  excitations.  This  can 
possibly  be  explained  by  noting  that  only  one  arch  support  is  vibrating 
under  unsymmetrical  excitations  which  in  turn  perturbs  any  sensitivities 
present  in  the  arch  one-half  the  amount  that  they  would  be  perturbed 
under  symmetrical  excitations. 

By  comparing  results  for  arches  with  various  a  and  R/H  values. 
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it  can  be  concluded  that  the  degree  of  agreement  between  theory  and 
experiment  is  governed  predominantly  by  the  value  of  a  and  to  a  much 
lesser  extent  by  the  value  of  R/H . 

It  was  calculated  theoretically  and  shown  experimentally 
that  under  symmetrical  excitations,  it  is  impossible  to  excite  unsym- 
metrical  modes.  However,  under  unsymmetrical  excitations  it  is  possible 
to  excite  symmetrical  modes.  This  can  be  explained  by  noting  that  un¬ 
symmetrical  displacements  are  obtained  by  superimposing  a  symmetrical 
displacement  plus  a  rotation. 

Although  out-of -plane  modes  were  not  studied,  they  were  noted 
to  occur  at  frequencies  having  values  approximately  16  times  the  fre¬ 
quency  of  the  first  flexural  mode.  This  agrees  with  the  observations 
of  Vol terra  and  Morel  1  (15). 


•v 


84 


Chapter  IV 
Closure 

4. 1  Summary 

This  investigation  was  concerned  with  the  development  of  a 
theory  which  could  be  used  to  investigate  the  dynamic  response  of  a 
thin  elastic  arch  subjected  to  cyclic  support  movement.  Only  flexural 
deflections  in  the  plane  of  the  initial  curvature  of  the  arch  were  in¬ 
vestigated.  The  arch  was  assumed  to  be  circular  having  a  constant 
mass  and  cross-section  and  the  center  line  was  assumed  to  be  inexten- 
sible.  A  closed  form  solution  was  obtained  for  the  homogeneous  part 
of  the  differential  equation.  The  steady  state  solution  was  assumed 
as  a  series  in  terms  of  the  free  modes  of  vibration.  Since  the  natural 
frequencies  were  spaced  sufficiently  far  apart,  it  was  found  that  an 
eight- term  series  gave  accurate  results  up  to  three  significant  figures 
for  arches  excited  near  the  fourth  mode.  Better  accuracy  was  obtain¬ 
ed  for  excitations  below  the  fourth  mode.  The  theory  presented  is 
applicable  to  arches  having  a  number  of  different  boundary  conditions 
such  as  pinned-end,  fixed-end,  free-end  or  any  combination  of  these. 

Experimental  evaluation  of  the  theory  was  presented  on  the 
basis  of  pinned-end  boundary  conditions.  Results  were  presented  for 
four  arches  having  values  of  121 .2  £  R/H  £  179.4  and  97.2  £  a  £  119.2 
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degrees.  Two  thicknesses  of  0.032  and  0.042  inches  and  a  width  of  one 
inch  were  used.  Other  arches  were  tested  in  the  above  range  of  R/H 
and  a  obtaining  equally  good  results,  however,  the  results  were  not 
included  in  this  report  as  the  four  arches  presented  cover  the  range 
of  R/H  and  a  tested  sufficiently. 

During  the  experimental  stage  of  the  investigation,  it  was 
discovered  that  arches  had  to  be  circular  within  one-half  their  thick¬ 
ness  before  reliable  test  results  were  obtained.  For  arches  not  per¬ 
fectly  circular,  the  experimental  amplitudes  were  much  greater  than  the 
theoretical,  indicating  that  some  parameter  was  highly  sensitive  in 
arches  having  a  large  variation  in  radius.  As  a  result,  the  arches 
did  not  respond  according  to  theory  where  the  assumption  of  circular 
arches  was  used.  Due  to  the  capacity  of  the  shaker  table,  arches  with 
a  <  90  degrees  were  too  stiff  and  could  not  be  tested.  It  was  found 
that  arches  with  a  >  125  degrees  were  too  susceptible  to  any  perturba¬ 
tions  and  were  not  tested  for  this  range  either.  On  the  basis  of  the 
series  of  tests  conducted,  it  was  noted  that  the  value  of  one-half  the 
opening  angle  predominantly  governed  the  degree  of  agreement  between 
theory  and  experiment.  The  agreement  was  governed  to  a  much  lesser 
extent  by  the  value  of  the  radius  to  thickness  ratio.  On  some  arches 
tested,  the  first  out-of-plane  mode  was  noted  to  occur  at  a  frequency 
approximately  sixteen  times  the  frequency  of  the  first  flexural  mode 
which  agrees  with  the  observations  made  by  V ol terra  and  Morel  1  (15). 


86 


4.2  Problems  for  Further  Study 

As  a  result  of  this  investigation,  other  problems  requiring 
further  study  are: 

1.  Determine  what  effects  concentrated  masses  attached  to  the 
arch  have  on  the  magnification  ratio. 

2.  Extend  the  present  theory  to  include  other  arch  shapes  such 
as  a  parabola,  cycloid,  or  a  catenary. 

3.  Investigate  the  dynamic  response  of  arches  for  large  radial 
displacements  i.e.  R+u  f  R. 

4.  Perform  experiments  to  test  the  inextensional  theory. 

5.  Investigate  the  out-of-plane  vibrations  caused  by  parametric 
excitations  (i.e.  when  plane  of  initial  curvature  of  the  arch  and  the 
plane  of  the  support  displacements  do  not  coincide,  parametric  excitation 
will  occur). 

6.  Investigate  further,  what  effect  variations  in  the  radius 
has  on  the  magnification  ratio. 
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